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JET DIFFUSION LIQUID GREATER DENSITY 


SYNOPSIS 


study the influence the density difference between jetand the sur- 
rounding fluid the diffusion circular vertical submerged turbulent jet 
presented. The results are summarized Figs. 1(a) and 1(b), which show 
how the density and the velocity along the axis the jet depend the initial 
density difference between the jet and the surrounding fluidand the velocity 
the jet when issues from nozzle. The theoretical basis for the families 
curves Fig. described together with laboratory experiments, which 
confirmed the theoretical considerations. The experiments were performed 
with water jets. The density the jet was 1,000 per (fresh water), 
the density the surrounding fluid (salt solution) varied from 1,020 1,050 
per Because the range densities covered the experiments the 
results apply the case sewage disposal the marine environment. 


INTRODUCTION 


For the case turbulent jet surrounded homogeneous still water, the 
degree dilution depends the quantity and velocity flow, the diameter 
the outlet-pipe, and the density difference between the fluid the jet and the 
surrounding water. The relationship between these factors described here 
for the case turbulent jet issuing vertically upward surrounding water 
with specific weight greater than that the jet. Theoretical information 


Note.—Discussion open until November 1960. extend the one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. June, 1960. 

Hydr. Engrg. Lab., Wave Research Projects, Univ. California, Berkeley, Calif. 


ork 
| 
q 
ne 
4 


presented together with results laboratory tests. The range density dif- 
ferences covered the experiments such that the results obtained apply 
the case sewage disposal the marine environment. 


THEORY 


the case homogeneous surrounding water withno currents, the velocity 


and the concentration point the jet dependin general the follow- 
ing variables: 


the distance from the nozzle the outlet-pipe measured along the 
axis the jet; 


the distance from the axis the jet; 

the diameter the outlet-pipe; 


the density difference between the fluid the jetat and 
the surrounding fluid; and 


The concentration defined as: 


which the density any point the jet and the density the sur- 
rounding fluid. 


Dimensional analysis that the velocity and the concentration de- 
pend the other variables as: 


and 
u 
and Reynolds number, 


The Reynolds number gives the ratio the inertia forces the viscous 
forces. Since this paper restricted turbulent jets, the Reynolds number 


considered sufficiently large have influence the functions 
Eqs. and 


and Brooks, Proceedings, ASCE, Vol. 86, No. March, 1960. 
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The Froude number gives the ratio the inertiaforces the gravitational 
forces. any point the jet the local basedon the local ve- 
locity and the local density-difference with the surrounding water large, the 
influence the gravitational forces may neglected that point. evi- 
dent that the inertia forces are less important greater distances from the 
nozzle. Hence possible that near the nozzle, the Froude number suf- 
ficiently large neglect the influence the buoyancy while may that 
greater distances from the nozzle the buoyancy has considered. 


According the problem may divided into the following 
parts: 


The non-buoyancy case, defined For this case the buoyancy 
negligible any point the jet. 

The buoyancy case, defined For this case the initial momen- 
tum the jet negligible. 

The intermediate case, defined finite value For this case the 
initial momentum predominant near the nozzle, the buoyancy larger dis- 
tances from the nozzle. 


discussion the various parts the problem follows. 

The Non-Buoyancy Case, regarding the non-buoyancy case 
Forstall and Hinze and van der Hegge Zynen,® 
Reid.8 

These studies yield the following results: 


For the zone flow establishment 


“On the Diffusion Heated Jets,” F.H. Schmidt, Sartoyck Tellus, No. 
August, 1957, pp. 378-383. 

Rouse, Transactions, ASCE, Vol. 115, 1950, pp. 639-664. 

“Momentum and Mass Transfer Submerged Water Jet,” Forstall and 
Gaylord, Journal Applied Mechanics, Vol. 22, No. June, 1955, pp. 161-164. 

“Transfer Heat and Matter the Turbulent Mixing Zone Axially Symme- 
trical Jet,” Hinze and van der Hegge Zynen, Seventh Internatl. Congress 
for Applied Mechanics, London, 1948, pp. 286-299. 

Experiments the Flow and Heat Transfer Heated Turbulent Air 
Jet,” Corrsin and Uberoi, N.A.C.A. Tech. Note 1865, April, 1949. 


Reid, Batelle Memorial Inst., The Rand Corp., Santa Monica, Calif., February, 


if 


= 
; 
Uo 2 
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For the zone established flow 


for 


C2; 


which 


the velocity the axis the jet c,, the con- 


centration the axis the jet, and Cg, and are dimensionless 
constants. obtain smooth transition between the zones that 


The experimental values for the constants and and the correspond- 
ing values and are shown Table 


The experiments Corrsin and Uberoi? show that the range 


constant. consequence the constants and not vary much for 


that is, within the range values that are interest sew- 

age disposal problems the ocean. 


d 
and 
and 
4 1 1 ‘4 


1.025 may describe the experimental results aver- 


age values for and according theory average values for and 
Since water jets are involved the disposal sewage, the following have 

These are the values obtained the Forstall and 


TABLE PUBLISHED DATA FOR NON-BUOYANCY CASE 


Type 

jet 

jets 0.080 

Forstall and water jets 

salt solutions 1.01 0.078 


Gaylord 


Keagy al. 


Hinze and 
van der Hegge 
Zynen 


mixture 
gasses 


Corrsin and 
Uberoi 


The Buoyancy Case the buoyancy case were made 
Morton, Taylor and Turner? and Rouse, C.S. Yih and 


The theoretical work Morton, al.,9 yielded the following results: 


and 


These formulas have been derived assuming that the velocity and buoyancy 
force are constant across the jet and zero outside and assuming that the 
jet had initial momentum. 


Gravitational Convection from Maintained and Instantaneous Sources,” 


Morton, Taylor and Turner, Proceedings, Royal Soc. London, Series 
March 1956, 11, pp. 1-23. 


“Gravitational Convection from Boundary Source,” Rouse, C.S. Yih and 
Humphreys, Sartoyck Tellus, Vol. No. August, 1952, pp. 201-210. 


The formulas are only valid for jet satisfying For jet with finite 
diameter the factor and must bereplaced the factor 2d). 

The constant must found experiment. Experiments made with water 
jets indicated that 0.093. 

possible write the equivalent Eqs. and for jet with finite 
diameter ina dimensionless form use the relationships Eqs. and 


-1/3 

and 


The experiments Morton, al.,9 showed that the distribution perpendicu- 
lar the axis the jet given the Gaussian curves: 


m e 


and 


Rouse, derived the 5/3 and 1/3 power Eqs. and theo- 
retically. Their semi-empirical formulas, based experiments done with air 
jets reduce dimensionless form to: 


5/3 
and 


The author’s experimental results agree better with Eqs. and 17, due 
Morton, al., than with Eqs. and 21, due Rouse, al., since the first 
set equations takes small initial momentum Therefore the 
author has preferred apply Eqs. and the following section. 

The Intermediate Case, Finite Values F.—The intermediate case has been 
studied theoretically Schmidt.3 The equations that derived suggest that 
the intermediate case may described approximately the equations for the 
non-buoyancy case (Eqs. 10, and the previously selected values 
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JET DIFFUSION 


and C4) from the nozzle certain distance from the 
nozzle, and the equations for the buoyancy case (Eqs. and 17) greater 
distances from the nozzle. The transition occurs the point where 
the curves that describe the non-buoyancy case and the buoyancy case inter- 
sect each other (Fig. 1). 

For sufficiently large values about 10) the ordinates the two 
sets points intersection satisfy 


with Cy, and according the previously selected values. 


For sufficiently large values and hence for sufficiently large values 


since for sufficiently large values the the zone flow es- 
tablishment the integral Eq. relatively small and then does not 


make much difference apply Eqs. and also for evalu- 


ate the integral. 

Eqs. and show, for sufficiently large values least, that the in- 
fluence the initial momentum larger than the influence the buoyancy 
force for while the contrary applies for This compatible with 
the description the intermediate case suggested the equations Schmidt 
which Fig. graphical representation. the following section experi- 
ments will described which the curves Fig. are confirmed. 

Fig. only applies for the case that the initial momentum and the buoyancy 
force have the same direction. For the case that the initial momentum and the 
buoyancy force have opposite direction, the distance approximation 
the distance from the nozzle, where the velocity along the axis will reduced 
zero. 


EXPERIMENTS 


The experiments check the families curves Fig. were made ina 
jet issued into the tank through pipes through the tank bottom with internal di- 
ameters 0.6 in., 0.8 in. and1.0 in. The end the vertical pipes was about 
in. above the bottom the tank. The upper edge the tank was constructed 
continuous weir. The discharge the jet was determined volumetric 
measurement the discharge over the weir. 

The density the fluid the jet was about 1,000 per The density 
the surrounding fluid varied from 1,020 per 1,050 per 
The surrounding fluid was solution crude table salt and water. 

Measurements the Density.—The density the axis the jet was mea- 
sured check the curves for the concentration. Samples the fluidof the jet 
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different points the jet with internal diameter 
3/16 in. The pipes were mounted frame that could moved through 
the fluid such way that the inlet opening described vertical line (Fig. 2). 
The discharge through the pipes was such that periodof sec was required 
collect sample 100 cc. The density the samples was measured with 
Westphal balance. The balance was accurate about 0.1 per 


TABLE 2.—RESULTS MODEL RUNS 
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The time history the density the fluidof the jet was measuredat differ 
ent points the jet. typical example shown Fig. Fig. 3(a) shows 
the time ordinary scale and shows the same data asemi- 
logarithmic scale. The first steep part the logarithmic plot (and the corre- 
sponding part the other plot) describes the establishment equilibrium. 


: 
d. in x Uo» in Cm Cc. ie 
0.6 23.8 +1/2 
10.6 40.7 1.6 3,360 195 2-1/2 
10.6 20.4 34.3 2.2 4,250 245 


The second less steep part describes how the equilibrium disturbed the 
fact that the surrounding fluid gradually mixed with the fluid the jet. The 
small dimensions the tank made impossible obtain real equilibrium. 
Therefore, the density corresponding the point intersection the two 
straight lines the logarithmic plot has been considered equilibrium 
density. Combined with any error made the density measurement itself, 
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(a) COMPARISON BETWEEN THEORY 
AND EXPERIMENTAL DATA 
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(b) DEVIATION BETWEEN EXPERIMENTAL DATA AND THEORY 


FIG, 4.—COMPARISON BETWEEN THE EXPERIMENTAL 
RESULTS AND THE THEORETICAL DATA 


this method will cause error about 0.5 per the determination 
the density difference between the fluid the jet and the surrounding fluid. 


Hence the possible absolute error the determination the dilution 
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The determining circumstances and the results the model runs regarding 
the density are shown Table For most the runs the Reynolds number 
sufficiently high assure turbulent flow conditions. 

The experimental results are plotted Fig. comparison between the 
experimental results and the theoretical data, based the theory Morton, 
shown Figs. 4(a) and (b). 

The experimental results agree well with theory. The approach, suggested 
the equations Schmidt, combine the results the non-buoyancy case 
with those the buoyancy case yields good description for the intermediate 
case. That Eqs. and 17, due Morton, al.,9 take small initial momen- 
tum account the reason for the goodagreement obtained for small values 

Measurements the Velocity.—A pitot tube has been used attempt 
measure the velocity-distribution the jet. The velocity readings were, how- 
ever, insufficiently accurate because the following reasons: (a) equilibrium 
conditions were not reached the tank during the model the velocity 
the jet function the distance from the nozzle, while pitot tube re- 
quires constant velocity between the static- and and (c) the 
small inaccuracies made the determination the velocity-head were im- 
portant because the magnitude the velocities involved (20-30 per sec). 

Since other instrument (small propellor, for example) was available, 
further attempt measure the velocities has been made. 


CONCLUSIONS 


The family curves shown Fig. 1(a) describes the dilution due sub- 
merged circular turbulent jet, issuing vertically upward higher 
specific gravity than the jet. The curves are based theory, which ade- 
quately confirmed model experiments. The density-range covered the 
experiments such that the curves may used when sewage disposal the 
marine environment considered. 

The family curves, shown Fig. 1(b), describes the velocity -distribution 
the jet. The curves are based theory. attempt has been made con- 
firm this set curves experimentally. 
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TRANSLATIONS FOREIGN LITERATURE HYDRAULICS 


Second Progress Report the Task Force List Translations the 
Committee Hydromechanics the Hydraulics Division 


INTRODUCTION 


attempt bringto hydraulic engineers up-to-date listing trans- 
lations foreign literature hydraulics, the intent force pre- 
report issue, fromtime totime, progress reports containing suit- 
able items. Interested readers are urged submit discussions (a) add 
this list, (b) offer suggestions for improvement, and (c) otherwise assist the 
task force infulfilling its aims. The present listis the second the 
publication, 1957, ASCE Manual Listof Translations For- 
eign Literature list appearedin the August 1959 Jour- 
nal the Hydraulics Division. This list considered addendum 


Manual and will includedin revision that Manual when the time 
deemed appropriate. 


ABRAMOVICH, “The Turbulent Jetin Moving Fluid.” Izv. Akad. Nauk, 
Otd. Tekh. Naul. 1957, No. pp. 93-101, 465. (On file the Dept. 
Scientific and Industrial Research, Lending Library Unit, London, England), 


AL’TSHULER, V., BAKANOVA, A., AND TRUNIN, F., “Phase Trans- 
formations During Compression Water Strong Shock Waves.” Akade- 


miya Nauk SSSR. Doklady, Vol. 121, No. pp. 67-69, 1958. Tr. No. 59-15750. 
(On file SLA). 


ANONYMOUS. “Symmetrical Cavitation Flow Arounda Flat Plate, Situated Be- 
tween Parallel Nauk SSSR. Otdeleniye Tekhnischeskikh 
Nauk, Izvestiya, No. 1946. Tr. No. 59-15095. (On file SLA). 


Note.—Discussion open until November 1960. Separate discussions should sub- 
mitted for the individual papers this symposium, extend date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. June, 1960. 
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‘Instrument for Measuring Cavitation During the Operation 
Water Turbines,” (Pribor dlyaizmereniya kavitatsii pri rabote gidroturbin). 
Izobretatel’stvo SSSR, Vol. No. 10, 40, 1957, DSIR 33866/CT. (On 
file Dept. Scientific and Industrial Research, London). 


ARKHIPOV, “The Formation Streaming Fluctuations Behind Solid 
Obstacle.” Dok Nauk SSSR, Vol. No. pp. 620-624, 1958. Tr. 
Soviet Physics Vol. III, No. published bythe Am. Inst. phys- 


DORFMAN, Rotating Rough Disc.” Journalof Technical Phys- 
ics the Academy Sciences the USSR, Am. Tr. Soviet Physics Tech- 
nical Physics, Vol. No. pp. 353-360, Publisher: American Institute 
Physics, Inc., Dec. 1956. 


DROBLENKOV, “Turbulent Boundary Layer aCoarse Curved Surface.” 
Izv. Akad. SSSR, Otd. Tekh, Nauk, pp. 17-21, 1955. 239. (On 
file the Dept. Scientific and Industrial Research, Lending Library Unit, 
London, England). 


DUBINSKII, “Vortex Izv. Akad, Nauk, Otd. Tekh. Nauk. No. 


pp. 3-10, 1955, C.T.S. 240. (On file Dept. Scientific and Industrial 
Research, Lending Library 


Flow Pastan Air Bubble aVertical Pipe.” Zeitschrift 
fur Angewandte Mathematik und Mechanik, Vol. 23, No. pp. 139-149, 1943. 
Tr. Olscher, 1958. (On file SLA: 59-10290). 


FEDOROV, motionof Plate Infinite Span Near the Free Surface 
Ideal Imponderable Fluid.” Reports Academy Sciences U.S.S.R. 
Vol. 116, No. pp. 373-376, 1957, 145. (Onfile Dept. Scientific and 
Industrail Research, Lending Library Unit, London, England). 


FERRARI, CARLO. “Wall Turbulence. Translations Section the mon- 
ograph Sulla Teoria della Turbolenza” (Turbulent Flow The- 
ory).” Fromlectures delivered the sessions Flow Theory” 
Varenna, Italy, Sep. 1-10, 1957. Tr. NASA, 1959. 


Circular Cross Section.” Journalof Technical Physics the Academy 
Sciences the USSR. Vol. 28, No. Tr. the American Institute Phys- 
ics, Vol. No. pp. 1681-1686. 


FRANKL’, “Energy Equations for the Motion Liquids Containing Sus- 
pended Solids,” Dokl. Adad. Nauk SSSR, Vol. 102, No. pp. 903-906, 1955 
C.T.S. 165. (On file the Dept. Scientific and Industrial Research, Lend- 
ing Library Unit, London, England). 


GOLIK, Z., and KARLIKOV, N., Connection Between Coefficient 
Viscosity and the Structure Substance the Fluid Akad. 
Nauk., Vol. 114, pp. 361-364, 1957, R.A.E. 711. (On file the Dept. Sci- 
entific and Industrial Research, Lending Library Unit, London, England). 


GRIGORASH, AND STAS’, “Fluid Surface Waves Caused Small 
Under-water Explosion Impulses. Determining Travel Time Tsunami 
Waves shore the Kurile Island. Remote Pressure Trans- 
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TRANSLATIONS 
mitter.” Published Academy Sciences (USSR) Marine Hydrophysical 


Institute (USSR), Vol. 11, pp. 18-27, 28-35, 42-55, 1957, Navy Tr. 1915. 


GUREVICH, Impact ona Lamina with Discontinuous Streamline 
Flow.” Prikladnaya Matematika Mekhanika (USSR), Vol. 16, No. pp. 116, 


1952. Tr. Gerta Cohen and Hirsh Cohen, 1959. (On file 
No. 59-15414). 


Fluid,” (Kolebaniya Plavayushchego Kontura Poverkhnosti Tyazheloy’ 
pages. Tr. Mildred Brode, David Fisher, and Alice 
Thrope, Feb. 1959. Order from: Library Congress (Photo dupl. Service) 
Literature Service Associates, Bound Brook, 


HERZOG, MAX. Dams” (Schalenformige Staumauern). Die Baut- 
echnik, Vol. 33, No. pp. 47-54, 1956. (Onfile SLA Translation Center, 
East Randolph St., Chicago 


IEVLEV, “Some Problems the Hydrodynamic Theory Heat Transfer 
the Flow Incompressible Liquid.” Dokl. Akad. Nauk, Vol. 86, No. 
pp. 1077-1080, 1952, R.A.E. 530. (On file the Dept. Scientific and In- 
dustrial Research, Lending Library Unit, London, England). 


JAROCKI, “Artificial Transverse Circulation” (Circulation transversale 
artificielle). Paper presented the Third Meeting the International As- 
sociation for Hydraulic Structures Research, pages, Grenoble, 1949. Re- 


port No. Tr. Margaret Petersen, March 1959. (On file WES 
Library). 


JUI, LI. “Water-power Resources and Conditions Hydro Development 
China.” Gidrotekhnicheskoye Stroitel’stvo (USSR), Vol. 27, No. pp. 20-26, 
1958. Order from: OTS, Dept. Commerce, Washington 25, 


KOZLOV, AND MOLOGIN, “On the Velocity Rise and Hydraulic 
Resistance Gas Bubbles Nauk SSSR. Otdeleniye 
Tekhnicheskikh Nauk. Izvestiya. No. pp. 1188-1197, 1951. Tr. Michael 
Roman, AEC, Tr. 3274, May 1958. (On file SLA). 


MEHAUTE, Perviousness Rock-fill Breakwaters Periodic 
Gravity Waves, Part (Permeabilite des digues enrochements aux ondes 
gravite periodiques, deuxieme partie). Houille Blance, Vol. 13, No. 
pp. 148-179, 1958, ATS-36L 31F, Associated Technical Services, Inc. 


MEHAUTE, “The Perviousness Rock-fill Breakwaters Periodic 
Gravity Waves, Part (Permeabilite des digues enrochements aux 
ondes gravite periodiques, Houille Blanche, Vol. 13, 
pp. 255-269, 1958, ATS-99L31F, Associated Technical Services, Inc. 


LEONTYEV, K., AND K., “Oscillatory Movements Shore- 
lines and the Formation Lagoons.” pages, April 1959. Tr. 
Hopper. ONI Translation No. 508, Office Naval Intelligence. 


LOGINOV, “Development and Progress Soviet Hydroelectric Power 


Construction over Years.” Tr. Gidrotekhnicheskoye Stroitel’stvo 
(USSR), No. 11, pp. 6-14, 1957. Order from OTS. 
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MAKSIMOV, “The Long-period Luni-solar Tide the Pacific Ocean.” 
Dok Nauk SSSR. Tr. in: Soviet Physics Vol. No. 1957, 
published the Am. Inst. Physics, 


MINSKII, Statistical Proof the Equations Filtration Motion.” 
Dok Nauk SSSR. Tr. in: Soviet Physics “Doklady,” Vol. No. June 
1957, published the Am. Inst. Physics, 


MINTS, Hydrodynamic Resistance Granular Layer Suspended 
Current.” Dokl. Adad. Nauk., Vol. 83, No. pp. 541-544, 1952, 834. 


(On file the Dept. Scientific and Industrial Research, Lending Library 
Unit, London, England). 


MISES, VON. “Computation Coefficients Outflow and Overfall-2) Efflux 
From Funnels.” 472 (partial translation) (Berechnung von Ausfluss und 
Ueberfallzahlen-2) Ausfluss aus Trichtern). Zeitschrift des Vereines deut- 
scher Ingenieure, Vol. 61, No. 22, June 1917. Tr. Keulegan, 
Bureau Standards, April 1959. (On file WES Library) 


NIKOLSKII, “The ‘second’ Form Motion Ideal Fluid Round 
Body” (Investigations Separated Vortex Flows). Dokl. Akad. Nauk, Vol. 
116, No. pp. 193-196, 1947, RAE 723. (On file the Dept. Scientific 
and Industrial Research, Lending Library Unit, London, England). 


OESTERHAUS, AND MULLER, ‘Plan for Fixing the Bed the Rhine 
between Strasbourg and Istein.” International Navigation Congress, 
Lisbon, September 1949. Section Inland Navigation, 2nd Communication. 
Tr. Margaret Petersen, Little Rock District, Corps Engineers, 
Aug. 1959. (On file WES Library). 


ized Hydrodynamics (Theory the Ultrasonic Interferometer).” Dok 
Nauk SSSR, Vol. CXXIII, No. pp. 821-825, 1958. Tr. in: Soviet Physics 
“Doklady” Vol. III, No. pub. the Am. Inst. Physics, 


ics.” Dok Nauk SSSR, Vol. CXVIII, No. pp. 468-472, 1958. Tr. in: 


Soviet Physics “Doklady,” Vol. No. pub. the Am. Inst. Physics, 


PANASENKOV, the Influence the Turbulence Liquid 
Jet its Atomization.” Zh, Tekh, Fiz. Vol. 21, No. pp. 160-166, 1951 
R.A.E. 409. (Onfile the Dept. Scientific and Industrial Research, Lend- 
ing Library Unit, London, England). 


PERELMAN, Forces Heterogeneous, Noncylindrical Pores, 
and the Problem Unidirectional Penetration Pores.” Dokl. Akad, Nauk., 
Vol. 85, No. pp. 153-156, 1952, C.T.S. 85. (On file the Dept. Scien- 
tific and Industrial Research, Lending Library Unit, London, England). 


POTAPOV, “New Method River Bed Hydraulics” (Novyy Metod Rus- 
lovoy Sochineniya, Vol. Chapter 1950. Tr. Nadia 


Georgevitch, Corps Engineers, Army, Washington, C., Mar. 
1959. 


REH, UWE. “Simplified Application Rehbock’s Formula for the Backwater 
Effect Caused Bridge Piers” (Vereinfachte Anwendung der Rehbock’schen 
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Bruckenstauformel). Die Wasserwirtschaft, Vol. 48, No. pp. 240-242, 
June 1958. Tr. Van Tienhoven, Army Engineer WES, Aug. 
1959. (On file WES Library). 


REINHOLD, TH., “The Mechanism Sedimentation Tidal Flats.” Mededel- 
ingen Geol. Sticht (N. S.), 1948-1949, No. pp. 75-80. Tr. Associated 
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RIEGELS, F., “The Flow Past Slender, Almost Axisymmetric Bodies.” Mitt. 
Max-Planck Inst., Stromungsforschung, Vol. 1952. Orderfrom: The Sec- 
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land. 


SELIVANOV, AND SHLYAPINTOKH, YA. “The Thermodynamic Pro- 
perties Air Thermal Ionization and the Shock Wave.” Zh. Fiz. Khim., 
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Horizontal Tubes.” Dokl. Akad, Nauk. SSSR. Vol. 104, Issue pp. 513- 
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search, Lending Library Unit, London, England). 
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FLOW LOSSES THE LOWER GILA RIVER 


Lawrence Pratt,! ASCE 


SYNOPSIS 


This paper describes method for estimating losses from infrequent flows 
the lower reaches the Gila River, Arizona. The same procedure could 


probably applied other intermittent streams sufficient data are avail- 
able. 


INTRODUCTION 


The usual procedure for studying analyzing flood flows stream 
wherein water flows continuously some kind mass-curve technique. The 
literature concerned with this problem voluminous. However, such proce- 
dure has little application intermittent streams because the available data 
are too discontinuous provide suitable universe for common statistical 
methods. 

Because the discontinuous nature the gaging records, the only ration- 
way study flood flows intermittent stream take them one 
time, not enmasse. gaging station intermittent stream, there will 
periods weeks, months, even years, when flow will recorded. Be- 
cause the many zeros any sortof average meaningless. Any 
given zero flow may reflect condition such that the last surface water may 
have disappeared only few feet, many miles, above the gage. there 
another gaging station upstream, obvious that the upper gage will record 


Note.—Discussion open until November 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. June, 
Cons, Engr., North Hollywood, Calif. 
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many more daily flows than the lower gage, provided that local inflow between 
the gages negligible absent. 

possible express the percentage loss between upstream and 
downstream station interms the quantity flow the upstream station and 
duration flow the downstream station. The procedure cumulate con- 
current flow quantities both stations, making proper allowance for time 
between the two. This allowance can most conveniently made 
matching the days peak flows, but not constant. When the streambed has 
been dry for long time, much water will required fill surface depres- 
sions and initiate infiltration. These requirements take time wellas wa- 
ter. When flow has been established, the time travel may reduced 
little third the time travel for dry streambed. 

flow continues, the percentage water lost betweenthe upper and lower 
stations changes, rapidly for the first few days and then more slowly, ap- 
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FIG, 1.—LOCATION MAP 


proaching limit. flow continues long enough indicate the value the 
limit, may inferred that this value primarily caused evaporation, 
Since equilibrium will have been approached between surface flow and infil- 
tration. Higher percentages loss represent the prior mentioned require- 
ments for filling depressions and initiating infiltration. 


GILA RIVER BASIN 


important intermittent stream southwestern United States the low- 
portion the Gila River, Arizona (see Fig. 1). The reaches considered 
herein are from Dome Gillespie Dam, since 1928, and from Dome Mc- 
Dowell and San Carlos, for 1904 and 1905. 

The Gila River Basin covers area 58,180 miles, which about 
88% Arizona, 10% New Mexico, and Sonora. The upper 
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reaches the Gila River and its principal tributaries are mountainous 
terrain within which are several peaks over 10,000 elevation. from 
this terrain that nearly all runoff derived. about 250 300 miles above 
the mouth the Gila River, several storage dams have been constructed 
impound this runoff. The larger dams are: Coolidge the Gila; Roosevelt, 
Horse Mesa, Mormon Flat, and Stewart Mountain the Salt; Horseshoe and 
Bartlett the Verde; and Carl Pleasant the Ague Fria. usable 
storage capacity reservoirs the Gila River and its tributaries was 3,446, 
000 acre-ft2. This storage capacity not quite twice the average annual in- 
flow the reservoirs. 

Downstream from the reservoir system the Salt River Valley contain- 
ing some 800,000 acres irrigated land and several cities and towns, the 
largest which Phoenix. nearly every year, all available surface wa- 
ter released from storage running unregulated tributaries wholly con- 
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FIG, 2.—FLOW GAGE FROM JULY 30, NOVEMBER 26, 1904 


sumed irrigation and domestic uses. addition, large quantity ground 
water pumped and consumed. The result that very little water appears 
below Gillespie Dam, the outlet the valley, and only rarely does any water 
reach the mouth the Gila River the gaging station Dome, miles 
above the mouth. The river distance from Gillespie Dam Dome about 
150 miles. Many the flows registered the Dome gage originate from in- 
frequent rains falling west Gillespie Dam return flows from irrigated 
areas between the two gages. 

Early brief period during the years 1904, and 1905, gag- 
ing stations were operated concurrently the Salt and Verde Rivers just 
above their junction McDowell, and the Gila River San Carlos (near 


Water the Gila River Basin and Adjacent Areas—A Summary,” open 
file report the USGS, 1952. 


USGS Water Supply and Irrigation Paper No. 133. 
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present Coolidge Dam) andat Dome, also called Gila City. During this period, 


the author, using data for other periods, estimate the average daily flows. 


there were storage reservoirs the basin. 

The Gila River was dry Gila City (Dome) from January 1904, (and for 
unknown time earlier) July 29, 1904. Flow the gage started July 30, 
and continued until November 26, after which date the river was again dry. 
Only gage heights are published for Dome, but rating curve was prepared 


Average daily flows are published for the three upstream stations. 

The apparent time travel flood peaks, after flow was established 
Dome, was seven days from McDowell Dome, and five six days from San 
Carlos Dome. While the flow passing Dome persisted for 119 days, 
time did become large. The maximum daily flow Dome was 4,560 day- 
sec-ft and the total flow was about 220,000 acre-ft. The combined flows pass- 
ing the three upstream stations was about 508,000 acre-feet. The difference, 
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FIG, 3.—FLOW FROM JANUARY JULY 19, 1905 


apparent loss, was 288,000 acre-ft. This difference ignors both unmeasur- 
inflows and consumptive uses, neither which can estimated with any 
assurance for this period. 

Flows McDowell and San Carlos were combined six days earlier 
than Dome. The combined flows were cumulated day day and the same cu- 
mulation was made for the flows Dome. The cumulative differences were 
considered represent the cumulative losses and were tabulated both day- 
sec-ft and percentages the combinedupstream flows. Percentage loss 
the end the flow period was 56.8%. Cumulative loss percentages between 
55% and 57% were consistent for more than the last half the period (see 
Fig. 2). concluded that this order loss percentage fairly represents 
the conditions existing 1904, with relatively low but persistent flow. 

The Gila River was dry the Dome gage from November 26, 1904, until 
January 16, 1905. Flow began that date and continued until July 19, 1905, 
except for February Very large flood peaks occurred this six 
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months period, reaching maxima 82,000 day-sec-ft February 95,000 
day-sec-ft March 20, and 64,000 day-sec-ft Except for re- 
latively small peak 13,000 day-sec-ft April 28, flows gradually decreas- 
from April 14, July 19, when the river was again dry Dome. Because 
the numerous large flood peaks, the calculated cumulative percentage 
loss not constant was that calculated for the 1904 period (see Fig. 3). 
For the last half the period the loss was fairly steady around 40%, but 
during the first half the loss varied from 21% 43%. The measured inflow 
for the period was about million acre-ft and the outflow Dome was about 
million acre-ft. Again, consumptive uses and unmeasured inflows were 
ignored. 

Later 1907, until recording gage was established 
May, 1929, there are records discharge Dome. recording gage had 
previously been established Gillespie Dam July, 1924. From May, 1929, 
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FIG, 4.—GILA RIVER, 1929 1955 


through Septernber, 1955, there were total 312 months. only those 
months were there recorded flows Dome excess 1,000 acre-fta month. 
For the same period there were months when the flow past Gillespie Dam 
exceeded 1,000 acre-ft. Therefore, two-thirds the months when there was 
appreciable flow past Gillespie Dam there was flow Dome. Under these 
circumstances, correlation monthly flows the two stations quite 
meaningless and correlation annual flows impossible. 

During the (1929-55) records studied the Gila River, there were 
only eleven occasions when flows passing Gillespie Dam reached the Dome 
gage for sustained periods seven days more. The longest periods sus- 
tained flow were 1941, for 114 days, and 1932, for days. Two flows, 
both 1931, continued for days and days. Seven other flows were sus- 
tained for days days. 


USGS annual water supply papers for the Colorado River Basin. 
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was found that the percentage the accumulated flow passing Gillespie 
Dam, but lost before reaching Dome, was related the duration flow 
Dome. was also found that the time travel was affected the antecedent 
condition the streambed. the bed had been dry for several months, the 
travel time was much eleven days. After the channel had been wetted, 
the time was reduced three four days. The time travel for each per- 
iod sub-period was determined matching the appearances flood peaks 
the two gaging stations. The differences between the cumulative flows 
Gillespie Dam and the cumulative flows Dome were found and expressed 
percentages the Gillespie flow. These differences were considered 
the cumulative losses the reach, ignoring any inflow caused local pre- 
cipitation. The determinationof any such local inflow wouldinvolve the rough- 
est kind estimates, there are records flow and insufficient num- 
ber and inadequate distribution rainfall records make useful rainfall- 
runoff correlations possible. Consumptive uses were small and were ignored. 
The results were plotted Fig. 

Fig. shows percentage loss drops quickly for the first few days 
and then flattens out, usually approaching limit, sometimes gradually and 
sometimes quite abruptly. For the two flows more than days duration 
the curves are hyperbolic form, which suggested the calculation em- 


pirical equationto represent the average these two occurrences. The equa- 
tion the form: 


which cumulative percentage loss, the first day percentage loss, 
days flow the lower gage, empirical exponent. For the 
reach from Gillespie Dam Dome the equation is: 


CONC LUSION 


seems probable that equation the same form could derived for 
any other intermittent stream sufficient data were available. The constants 
would undoubtedly different. 

For those streams wherein flows are more frequent and longer average 
duration than for the lower section the Gila River, relationship linking 
duration flows channel losses, flows the lower station, would 
have more significance. the economic value the water couldbe establish- 
would then possible determine what extent salvage should 
undertaken, storage channel control, increase the duration flows 
and thereby decrease the cumulative percentage loss. 


ACKNOW LEDGEMENT 


The Gila River studies were made the author while was employed 
the Colorado River Board California. Fig. with explanation its de- 
rivation, was presented the United States Supreme Court the case 
Arizona vs. California, al, Elder, ASCE, Hydrographic Engin- 
eer, the Metropolitan Water District Southern California. 


| 
s 4 
By 
b 
; 
tat 
\ 


2524 June, 1960 


Journal the 


HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


CAPACITANCE METHOD MEASURING WATER FILM 


Black! 


SYNOPSIS 


water film thickness gage was developed for making rapid measurements 
the depth water sheet flowing plane The gauge uses the 
electrical capacitance between fixed metal plate and the water surface 
sense the air gap and, use differences, measure the water depth. 

Rapid readings the depth disturbed water sheet flowing down plane 
were made, The dynamic range the instrument was 0.130 with accu- 
racy +0.003 in. Transverse roller waves in. long were resolved with this 
instrument; through the use special “probe,” waves 1/2 long were re- 


INTRODUCTION 


The Naval Radiological Defense Laboratory currently studying the 
transport radioactive particulate contaminants water Stream depth 
one the most important factors the equations describing sheet-flow 
water incline This depth can measured precisely and accu- 


Note.—Discussion open until November 1960, extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
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Chem. Engr., Aerojet-General Nucleonics San Ramon, Calif. Formerly with the 
Radiological Defense Lab., San Francisco, California. 

Burge, Naval Radiological Defense Lab. Review and Lecture, USNRDL-RL- 
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rately in. the use micrometer point gauge when stream flow 

uniform and the surface glassy smooth. However, when water flows down 

incline plane with slope greater than pronounced transverse roller 

waves are generated short the entrance the plane, and these 

waves introduce serious problems inthe use point gauge for measuring 

depths. Wave height not constant for given flow condition and trough 

measurements are questionable, the “end point” for determing contact 

subjective. Furthermore, information concerning wave-shape results from 

the use this techique. 

review the literature revealed two other methods determining water 

sheet thickness. These were optical which requires the water 

flowing the outside cylinder, and capacitance utilizing 

frequency-modulated The first method was impossible due the un- 

suitable geometry the apparatus required tobe usedand the secondwas con- 
sidered likely subject electronic drift. 

future progress the experimental study sheet flow depends partly 
the development instrument that can used for making depth meas- 
urements during non-uniform flow, the objective this project was develop 
instrument capable making depth measurement sheet water flow- 


ing over plane surface, and obtain adequate response demonstrate wave 
shape. 


APPROACH 


There were several potentially useful means for determining the thickness 
water sheet, grouped under two basic approaches, The first was utilize 
absorptive properties relating the thickness the water sheet (a) de- 
tecting the absorption beta radiation passing through the water, (b) de- 
tecting the absorption light passing through dyed water. The second basic 
approach was measure the the air-water interface the 
thickness (a) subtracting the distance the stream bed detecting con- 
surface electrically optically, (b) detecting the beta backscatter 
radioactive source fixed above the surface, (c) detecting the electrical 
capacitance metal plate fixed above water and the (electrically ground- 
ed) water surface. 

brief consideration these possible means from the viewpoint utiliz- 
ing readily available equipment revealed that most them had shortcomings. 
Methods dependent radioactive sources would require relatively expensive 
equipment and heavily shielded radiation Electrical optical-contact 
detection equipment was Light-absorption methods, utilizing dyed 
water, would necessitate recirculating system anda translucent trans- 
parent stream bed, and would complicate delay other studies requiring 
transparent water sheet. However,a was commercially 
available utilizing patented mechanic-electric transducer which had spec- 
ified sensitivity output per change capacitance and less than 

0.1% drift The detection electrical capacitance between fixed ref- 


Transfer Condensing Vapors,” Kirkbride, Transactions, Amer. 


Inst. Chem. Engrs., 30:170, 

gelin, Chem. Engr. Progress, Vol. 48, No. 11, 1952, 557. 
“Mechanical Electric Transducer,” Kurt Lion, Rev. Sci. Instruments, 27:222, 
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erence plate and the water surface would avoid the disadvantages contact 
with the water dyeing the water. 


APPARATUS 


The commercially available capacitance transducer essentially capac- 
itance comparator andgives adirect current (D. that proportional 
the percentage difference between the values two capacitors witha common 
The transducer excited and permanently connected 250 
carrier generator. One the capacitors, plate, the “probe,” fixed 
above and parallel the electrically grounded water. The other, C2, ad- 
justable reference capacitor. The relative positions these and the other elec- 
trical units are shown Fig. They are described the Appendix. The 
filter and power supply were constructed according specifications 
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the manufacturer the capacitance transducer. The amplifier and 
readout either anoscilloscope with acameraor recording oscillograph. 
The probe area was arbitrarily chosen 1.0 sqin. (6.42 
designed vary with water-sheet thickness changing 0.130 in. 
(0.330 cm), and have mean capacitance, including fixed stray capacitance, 
between the probe and the water surface considered the same that 


between two parallel plates first approximation setting the probe para- 
meters; 


which the capacitance denotes probe area (cm2), and 
the distance from probe water (cm) 
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With the values given previously, two equations can established: 


6.42 
6.42 
Cy + 0.3 = 0.0885 


Solving these equations gives for and 0.38 in. for 
The probe (which removeable) and shield, that encloses and fixes the 
stray capacitance the transducer, reference probe lead wires 
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FIG, 2.—APPARATUS 


are rigidly fastened the probe mount, This assembly has been offered re- 
cently unit the transducer manufacturer. The probe mount held ina 
vise which rides vertical calibrated micrometer screw. With this mi- 
crometer screw, the probe can raised lowered known distance from 
the stream bed. The instrument shown Fig. 
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WATER FILM 
CALIBRATION AND PERFORMANCE 


reasonable compromise between adequate sensitivity and acceptable 
dynamic range was probe 1-1/4 in. wide 3/4 in. long, designated 
the “standard” probe. With the set 0.400 in. above the stream 
bed, and the differential amplifier gain set give full-scale response 

1-v signal, the dynamic depth range 0.030 0.160 in. 

The maximum sensitivity the instrument for measuring change depth 
approximately 0.0003 in. This sensitivity can used for measuring water 
surface undulations; however, accuracy not great enough use the full sen- 
sitivity effectively determining water film thickness. The sensitivity the 
instrument utilizing the 0.130 in. dynamic range +0.001 in. 

The instrument was calibrated place with the water flowing. 
metal plate known thickness (and approximately the thickness the 
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FIG, 3.—TYPICAL RECORDINGS ILLUSTRATING 
CALIBRATION AND OPERATION STEPS 


water sheet) was placed under the probe displace the water and give mid- 
scale reference reading. The probe was lowered advancing the micrometer 
screw until the output meter barely read off-scale; then the probe was raised 
0.010 in. increments (or other convenient intervals) until the opposite scale 
limit was reached, and incremental readings were recorded provide scale. 
Following this, the probe was re-set the mid-scale reading, and the reading 
for the plate was recorded provide reference, thus, completing the cali- 
bration. The reference plate was removed, and the instrument was then meas- 
uring the water sheet thickness. Fig. shows specimen ofthecalibration 
data and depth data. 

The accuracy the instrument was determined comparing its readings 
with water depth measurements smoothly flowing water taken with mi- 
crometer point gauge. Agreement was 0.003 in. 
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The resolution the instrument longitudinal waves using the standard 
probe was measured replacing the water sheet with flat metal plate. This 
plate had several widths slots cut adepthof 0.030 in. represent “waves.” 
This plate was moved under the probe, place the water sheet, and the slot 
depth, indicated the instrument, was recorded. The results this test 
tude, and 1.5 slot isthe shortest that can resolved ampli- 


tude. This resolution assumed the same the instrument resolution 
longitudinal water waves. 
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Aspecial probe was constructed shorter waves. probe 
short 0.080 in. diameter rod fixed perpendicularly tothe stream 
bed and shielded along its length with 0.25 in. diameter tube. The results 
resolution tests for this smaller probe are also presented Fig. in. 
slot can resolved with 50% response amplitude. This probe has low 
sensitivity nonlinear response water-sheet thickness, and, therefore, 
was not generally used for thickness measurements. However, can cali- 
brated the same manner the standard probe. 


CONCLUSIONS 


The accuracy the instrument, determined comparison with the mi- 
crometer point gauge, limited the conductivity the water. Although the 
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WATER FILM 


impedance the transducer high relative that the water, slight fluc- 
tuations the conductivity the water due depth and ion content affect the 
resistance and, therefore, the apparent “location” the grounded side 
Simple tests indicated that accuracy canbe improved +0.001 in. increas- 
ing salinity approximately 0.01 Normal 

The resolution the instrument waves is, part, limited the dimen- 
sions the probe. The shortest logitudinal wave that can resolved equal 
the length the probe. 

When the small probe used, the resolution waves limited diver- 
gence the force from the smallarea the probeto the 
water surface. the electrostatic lines could made parallel applying 
the proper voltage electrostatic ring around the probe, then, the resolu- 
tion the instrument would limited only the smallest probe which could 
supply readable signal. 

The instrument sensitivity changes water depth limited the elec- 
trical noise level (which approximately millivolts), the maximum probe 
area wave resolution requirements, and the minimum air gap 
between the probe and the water surface, determined the anticipated fluc- 
tuation depth. 

Instrument stability excellent; drift was within +0.001 in. over 4-hour 
periodaftera Thetransducer life appears reasonable; 
the unit required repair after one year modest use. 

Instrument sensitivity, useful range, and wave definition are interrelated, 
indicated above, must balanced against each other for particular appli- 
manipulating the value the area and shape the probe, the 


distance the probe from the mean water surface, and the type readout 
recording equipment, this depth-measuring instrument can applied va- 
riety hydraulics studies involving surface. The rapid indication 
water surface and ease reading, even with disturbed surface, can save 
much time. Theavailability electrical signal output suggests the possibilities 
centralized simultaneous recording multiple stations andautomatic control. 


APPENDIX 


Equipment used for this report: 
Capacitance Transducer; Delta Unit Model 901-1, Decker aviation Corp., Phil- 
adelphia, sensitivity ufd initial capacity. 


Regulated Power Supply; NRDL design (vacuum tube regulation) comply with 
specifications furnished manufacturer capacitance transducer. 


Filter; NRDL-fabricated per diagram supplied capacitance transducer in- 
struction manual, 


Direct-writing Oscillograph; Sanborn, Model 128, 200 mv/cm sensitivity, dif- 
ferential input. 


Oscilloscope; Dual Beam, Type 322A, DuMont Lab., Clifton, 200 
mv/cm deflection sensitivity, differential input (one beam required). 


Oscillograph-record Camera; Type 297, DuMont Lab., Clifton, 
Constant Voltage Transformer; Solavolt No. 50106, 500 Watt. 
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SEDIMENTATION ASPECTS DIVERSION OLD RIVER 


Fred 


SYNOPSIS 


recent years the Old and Atchafalaya Rivers have diverted progres- 
sively increasing percentage Mississippi River flow. The damaging conse- 
quences posed this continuing action necessitated controlled diversion. This 
paper presents the sedimentation aspects channel evolution leading ne- 
cessity shown that the distributary system state im- 
balance sediment transport capacity, the nature which led conclusion 
that the flow capacity the Atchafalaya River can expected annu- 
ally regardless the type water year. This conclusion substantiated 
field observations the Simmesport discharge range. 


INTRODUCTION 


The Atchafalaya River, which connected with the Mississippi River Old 
River, the farthest upstream distributary the Mississippi River. Fig. 
plan map showing the Mississippi River below latitude Natchez, the lower 
end the Black and Red Rivers, the Atchafalaya River, and connecting Old 
River. From the head Old River approximately 140 miles the Gulf 
Mexico way the Atchafalaya River. way the Mississippi River 
past Baton Rouge and New Orleans 320 miles, more than twice the dis- 
tance via the Atchafalaya River. Old River approximately 6.5 miles length. 


Note.—Discussion open until November 1960. extend the closing date one month, 
request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. June, 1960. 

Chf., Hydr. Model and Sedimentation Sect., Branch, Army Engr. Div., 
Lower Missippi Valley, Vicksburg, Miss. 
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The lower miles the Red River channel and the lower portion Black 
River, whose principal flow from the River, are backwater 
basin containing numerous shallow lakes and low-level swamp areas. the 
peak Mississippi River floods, and coincident peak flows the Atchafalaya 
River, the water-surface slopes through this backwater area are very flat. 
Consequently, the Red and Ouachita Rivers system can contribute only very 
fine sediments the Atchafalaya River high stages. 

recent years the Old and Atchafalaya Rivers have diverted progressively 
increasing percentages Mississippi River flow. Studies this development 
led the conclusion that this diversion from the Mississippi River will reach 
“critical” stage (approximately 40% the major river’s flow) 1975. The 
consequences such eventuality necessitated controlled diversion. Con- 
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STUDY 


FIG, 1.—RIVER SYSTEM BELOW LATITUDE NATCHEZ 


struction control structures was initiated 1955 and now essentially 
complete. All work, including new channels, bank protection, lock, and the 
closure Old River, scheduled for completion the end 1963. Details 
this project which include river history, the general problem, hydraulic re- 
quirements control works, type structures and foundation design, have been 
presented previously.2 

Previous discussions this problem area have dealt mostly with flow and 
channel dimensional changes. This paper presents the sedimentation aspects 
channel evolution leading tothe necessity controlled diversion. The sedi- 
ment transport characteristics channel reaches the provide 


River Diversion Control: Symposium,” Transactions, ASCE, Vol. 123, 1958, 
pp. 1129-1181. 
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SEDIMENTATION 


excellent example sediment transport capacity imbalance interconnected 
system rivers. For the purposes of:this paper, the discussion will con- 
fined Old River and the upper miles the Atchafalaya River. 


ATCHAFALAYA RIVER 


The Atchafalaya River flows, for the most part, through the leveed Atcha- 
falaya Basin Floodway which averages about miles width. project de- 
sign flood stage this floodway will required pass 1,500,000 cfs which 
50% the project design flood flow this latitude. The river channel 
the approximate center the floodway and leveed from Mile 53, meas- 
uring from its head. This leveed from 2,500 3,500 width 
and flanked the west side the West Atchafalaya Floodway and the 
east side the Morganza Floodway. From Mile 115, flow channeled 
through low-lying swamp and shallow lake areas. 

virtue lesser distance sea level, complemented widespread 
flow the lower basin, the water surface slope the upper leveed reach 
appreciable flow steeper thanthat the Mississippi River below Old River. 
Also, the leveed reach rather deep comparison with its width. For in- 
stance, 1957 the bankfull stage width-to-depth ratio the Atchafalaya River 
Simmesport was compared the Mississippi River just below 
the head Old River. Old River Torras, 1.1 miles from the Mississippi 
River, itwas 39. Mean depths measured bankfull stage were approxi- 
mately 60, 50, and ft, respectively, Simmesport, Torras, and Red River 
Landing. 

cannot stated being exact uniform occurrence, but there 
evidence that degradation occurs inthe Atchafalaya River anytime medi- 
and higher stages coincident with low velocities Old River. these 
stages Red River sediment supply the sand range negligible because 
flat slopes the backwater area. Under any circumstance the sand fraction 
Red River sediment much finer than that coming from the Mississippi 
through Old River. Therefore, the important sediment supply which affects 
the regimen the Atchafalaya River comes from Old River. Conversely, when 
there appreciable flow the Atchafalaya coming almost entirely from the 
Mississippi River, some aggradation the Atchafalaya expected. 

Extreme cases the latter condition are infrequent, the last occurrence 
being June, 1943. that time flow the Atchafalaya was excess 
400,000 cfs and all was coming from the Mississippi via Old River. Av- 
erage velocities River that time reachedas high fps compared 
with less than fps the Atchafalaya. Roughly, this would make the supply 
sediment the Atchafalaya about 50% excess its transport capacity. 
Observations indicate that bankful stage the mean bottom elevation the 
Simmesport discharge range may raise lower from ft, depending 
the percentage Atchafalaya River discharge being contributed through 
Old River. 

1952 study was made bed material transport capacity the Atcha- 
falaya River from its head Krotz Springs (Mile 41). This study was based 
1949-50 channel conditions and water stages observed during the 1950 flood. 
These data showed that water surface slopes this reach were not uniform 
either with respect stage location, nor was the channel section uniform 
average width and depth. the downstream direction the channel decreased 


e 
i 
sar 


rg 


width, increased depth and, general, there was increase water 
surface slope. was indicated that the slope also increased with increasing 
stage. 

study the effects which these varying channel characteristics might have 
sand transport, the river was divided into three reaches shown Fig. 
This division reaches, Barbre Landing Simmesport, Simmesport Mel- 
ville, and Melville Krotz Springs, was dictated location gaging stations 
which daily observations are made. 

Sand transport rating curves were then developed for each these reaches 
general accordance with the bed-load function developed 
that time the bed-load function theory was relatively new and 
initial studies were directed towards determining its applicability the Atcha- 
falaya River. Sediment data from frequent observations the Simmesport 
discharge range were available for this purpose. was found that application 
the principles the bed-load function theory, with afew modifications, would 
produce results reasonable agreement with the measured sediment loads. 

shows the sand transport rating curves developed this method for 
the three reaches. interesting note that the rating curves show the 
leveed reach channel be, general, state equilibrium flow 
400,000 cfs. This approximately bankfull flow the leveed reach. Above 
400,000 cfs, the transport capacity increases the downstream direction. The 
difference between Reaches and more less constant, but Reach 
itis increasingly greater increases. Thus, indicated that 
that time there was sediment transport capacity imbalance the upper 
reaches the Atchafalaya River stages above bankfull. 


OLD RIVER 


Old River approximately 6.5 miles length and the connecting channel 
between the Atchafalaya and Mississippi Rivers. Its mean depth bankfull 
Mississippi. Formerly its width was several hundred feet less than the Atcha- 
falaya, but recent years its rate enlargement has exceeded that the 
Atchafalaya and now there little difference between the two, this respect. 
The amount sand supplied the Atchafalaya River from outside sources 
includes that which passes into Old River from the Mississippi, plus and minus 
fillor scour River, plus that contributed the Red and Ouachita Rivers 
system. The latter system not considered important contributor ma- 
terial that will affect the regimen the Atchafalaya River. 

Discharge through Old River not function stage; itis, rather, 
dependent relative stages its head and mouth. fast rising Mississippi 
coincident with small contribution from the Redand Ouachita system produces 
high velocities River. These velocities are sometimes much higher than 
those attained the Atchafalaya the same time. Conversely, rapidly fall- 
ing Mississippi following flood crest results low velocities Old River 
with sustained high velocities the Atchafalaya. Old River flowis also affect- 
appreciable flow from the Red and system. 1945a large con- 
tribution from the latter system coincident with falling Mississippi resulted 
reverse flow Old River. This was not uncommon prior 1945. 


Dept. Agric., SCS, Tech. Bulletin 1026, September, 1950. 
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FIG, 2.—SEDIMENT STUDY REACHES 
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few instances since that has been reduced tovery small amounts, but 
there has been recurrence reverse flow. This attributable the in- 
creased flow capacity the Atchafalaya River. 


SEDIMENT TRANSPORT 


noted previously, the interrelated river reaches under discussion pro- 
vide example sediment transport capacity imbalance. present 
the magnitude imbalance necessary state sediment loads quanti- 
tatively. Also, order show the degree imbalance for full range 
water-type years, number years must considered. This involves la- 
borious computations carried out great detail. 

was not considered necessary for purposes illustrating this imbalance 
include the entire bed material load, nor was considered necessary 
make detailed load computations. Mechanical analyses number bed- 
material samples taken the study reach showed medium sand (0.25 0.50 
mm) the dominant size, comprising approximately 35% the total. Since 
any grain size group would serve show relative transport capacities quanti- 
tatively, determinations were made only for the medium sand fraction. For 
analytical load determinations, was found that close approximation the 
results detailed computations could attained simple formulas requiring 
only knowledge the average velocity depth. These formulas, applied 
velocities and depths, measured frequent intervals the Torras and 
Simmesport discharge ranges, are the bases for the quantitative medium sand- 
loads presented. 

examination mean annual discharges the principal contributary 
rivers led the adoption three water-type years for study. These were 
classified according tothe magnitude average annual flowin the Mississippi 
River high-, medium-, and low-water years. For each these types three 
years were selected representative various combinations mean annual 
flows the contributary streams. The years 1950, 1945, and 1951 were se- 
lectedas representative high-water years; 1957, 1958, and 1946as medium- 
water years; and 1944, 1953, and 1954 low-water years. 

the high-water years the average daily Mississippi discharge was from 
826,000 cfs (1951) 879,000 cfs (1950); for medium-water years from 605,000 
cfs (1946); 685,000 cfs (1957); and for low-water years from 333,000 cfs 
(1954) 567,000 cfs (1944). Variation average daily contributary 
the Red and Rivers system, that has large influence flow Old 
River, was from 46,000 cfs (1951) 102,000 cfs (1945) for the high-water 
years; 70,000 cfs (1958) 88,000 cfs (1957) for the medium-water years; and 
from 25,000 cfs (1954) 55,000 cfs (1944 and 1953) for the low-water years. 

Although 1944 was placed the low-water year category, borderline 
case from the standpoint Mississippi and could just well have 
been placed the medium-year category. However, from the standpoint 
and Atchafalaya Rivers flow, comparable with 1953, which 
true low-water year. tabulation average daily discharges and com- 
puted sediment transport capacities for the years selected for study shown 
Table 

Since sediment transport capacity directly related stream velocity, 
comparison daily average velocities Simmesport and Torras for the se- 
lected years shown Figs. Fig. shows the high-water years 


ts 
| 
aye 
q 
i 
: 
¥ 
j 
4 
va 


SEDIMENTATION 


1950, 1945, and 1951. each these years the annual medium sand transport 
capacity the Atchafalaya River exceeded that Old River. The greatest 
difference was 1945 and amounted nearly million tons. The next highest 
was 1950 with nearly million tons, followed 1951 with one million tons. 
comparison the velocity graphs shows why the difference was much 
greater 1945 than 1951. 

The volume water passing Natchez 1951 was about equal that 
1945, but the peak discharge was much higher This reflected the 


TABLE DAILY DISCHARGES AND COMPUTED SEDIMENT 
TRANSPORT CAPACITIES 


Average Daily Discharge Medium Sand? Excess Ca- 
pacity River over 
Million Tons Old River 
Missis- Atchafa-| (Medium sand 
sippi Atchafa- laya fraction 


Period Riverb 


laya Tons) 


2,832,000 
3,837,000 
1,030,000 


1957 531,000 
Apr-3 Aug 4.058(4) 2,824,000 
Oct-31 Dec 1.168 705,000 
1958 6.192 2,098,000 
1946 2,238,000 


1944 


1,893,000 


1953 
May-31 Jul 1,565,000 
Aug-31 Dec -61,000 
1954 


242,000 


Bed material samples show the average medium sand fraction about 35%. 
This percentage was used determination transport 
Above Old River. 


tributaries the backwater area. 


hydrographs shown Fig. for the Atchafalaya River, the 1945 peak being 
630,000 cfs and 1951, 425,000 cfs. Also 1945 there were sustained high 
flows from the system which was deterrent flow through 
Old River. Consequently, Atchafalaya River velocities were greater for 90% 
the period from January August 17th. Old River velocities were domi- 
nant for the most partonly during the low-water season which time sediment 
loads are relatively small. 

The year characterized sustained periods near-bankfull flow 
the Mississippi River. the months its flow equalled exceeded 


A 
1950 879 297 195 5.699 8.531 
1945 832 264 147 102 3.396 7.233 
1951 826 256 195 3.925 
Medium Water Years 
igh” 
Low Water Years 
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FIG, 5.—VELOCITY RELATION—MEDIUM WATER YEARS 
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SEDIMENTATION 


1,000,000 cfs. Flows less than 300,000 cfs occurred only days during the 
year. contrast, the year 1945, there were only months which the flow 
equalled exceeded 1,000,000 cfs, and there were days when the flow was 
less than 300,000 cfs. the year 1951 there were sustained periods ap- 
preciable the Red and system. Therefore, flow Old River 
was generally steady, and its mean stream velocity compared favorable with 
that Simmesport throughout most the year. result, the excess sedi- 
ment transport capacity the Atchafalaya River was much less than 1945. 

The year 1950 was characterized more than usual dominance Old River 
velocity during both the spring rise and the low-water season. However, this 
was more than offset greater velocities the Atchafalaya River during most 
the high-water period. The superior velocities and large volume flow 
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FIG, 6.—VELOCITY RELATION—LOW WATER YEARS 


the Atchafalaya River this time accounted for over 50% its 2,832,000 tons 
total excess for the year. 

Fig. shows comparative velocities the medium-water years 1957, 1958, 
and 1946. each these years the annual medium sand transport capacity 
the Atchafalaya River exceeded that Old River. 1957 the higher ve- 
locities the Atchafalaya River May, June, and July were sufficient per- 
mit the Atchafalaya River sediment-transport-capacity during this period 
alone exceed that for the entire year Old River. This illustrates the po- 
tential circumstances for great imbalance sediment transport capacity. 
This year also characterized the abnormal situation higher velocities 
the Atchafalaya River during most the low-water period. both 1958 and 
1946 the excess capacity the Atchafalaya River was more 
the medium sand fraction. 


. 
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Fig. shows comparative velocities the low-water years. other 
water-type years, the Atchafalaya River’s capacity for transporting sand ex- 
ceeds that Old River each these cases. Here again the potential for 
great imbalance illustrated. For example, May, June, and July, 1953, 
the transport capacity the Atchafalaya River was nearly times that Old 
River. Further, those months the transport capacity the Atchafalaya 
was over times that Old months, August through December, 
low-water period during which velocities Old River consistently exceeded 
those the Atchafalaya River. 

1944 the excess transport capacity the Atchafalaya was comparable 
some medium-water years, and nearly twice that high-water year 1951. 
Year 1954 was exceptionally low water year. Even so, and spite su- 
perior velocities Old River during most the year, the Atchafalaya River 
showed slightly greater annual capacity for transporting medium sand. The 
balance infavor the Atchafalaya River for this year was effected the May 
Mississippi River rise concurrent with increased flow from the Red and 
system. From the first May the fourth June the sediment 
transport capacity the Atchafalaya River exceeded that Old River ap- 
proximately one million tons. From June October the Old River ca- 
pacity exceeded that the Atchafalaya approximately 300,000 tons. This 
shows, again how effectively short period higher velocity and discharge 
the Atchafalaya River can offset long periods lower discharge which are 
favorable superiority transport Old River. 


CONCLUSIONS 


conclusive from the data presented that the distributary system the 
Mississippi River the vicinity the latitude Old River state ofim- 
balance sediment transport capacity. This imbalance has exist 
within the leveed portion the Atchafalaya River itself during large floods, 
and more importantly its relation with supplying streams all types 
water years from low high, when considered annual basis. 

The end result this state imbalance can lead only one direction— 
continuous annual enlargement the Atchafalaya River regardless water- 
type year. substantiation these findings from the standpoint sediment 
transport capacities, matter record that the rating curve for the 
Simmesport discharge range, during the period 1951-58, lowered ap- 
proximate bankfull discharge. this seven-year period there were high- 


water years; three the years were medium-year magnitude and four were 
low-water years. 
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DRAG AND LIFT SPHERES WITHIN CYLINDRICAL TUBES 


Donald Young! 


SYNOPSIS 


method for evaluating the lift and drag forces that are exerted mov- 
ing fluid spherical particle resting the bottom cylindrical tube 
described. Experimental data are presented show the variation the lift 


and drag force with the particle pipe-diameter ratio and the Reynolds num- 
ber. 


INTRODUCTION 


One the most important aspects the general problem the flow 
suspensions the study the conditions necessary for maintenance the 
solid particles suspension. This specific problem is, turn, related 
the evaluation the fluid forces that act the individual particles. The 
sultant dynamic fluid force acting particle can resolved into two 
ponents, drag force parallel the direction flow, and lift force normal 
the direction flow. 

Although some attention has been given the effect the drag component 
the behavior suspensions, little work has beendone connection with the 
determination the lift component. probably true that certain prob- 
lems the lift force has little significance, however, for the study incipient 
motion particles resting the bottom channels pipes the lift force 


Note.—Discussion open until November 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. June, 1960. 

Assoc. Prof., Dept. Theoretical and Applied Mechanics, Iowa State Univ. 
Science and Tech., Ames, Iowa. 
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may primary importance. Leliavsky has summarized2 the work that 
has beendone inthis area and points out that “an almost unexplored aspect 
the problem.” The purpose the investigation reported herein was 
develop method for evaluating the lift and drag forces that act single 
spherical particle located the bottom tube, and study the parameters 

that affect the magnitude these forces. 


MEASUREMENT LIFT AND DRAG FORCE 


The free-body diagram for spherical particle contact with fixed 
force has been resolved into lift component and drag component 
should noted that, general the resultant fluid force does not pass through 
the center gravity the particle. White has investigated3 the loca- 
tion the drag component some detail. The weight the particle minus 


FIG, 1.—FREE-BODY DIAGRAM 
the buoyant force isdenoted and theforces Fand com- 
ponents the boundary reaction the particle. 
Equilibrium equations for the particle are 
“An Introduction Fluivial Hydraulics,” Leliavsky, Constable and Co., 1955, 
64. 


Royal Soc. London, A174, 1940, 
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for moments with respect axis through the center gravity the par- 
ticle diameter 

varying the rate flow past the particle, the drag component can 
changed and the particle made move down the plane with constant 
velocity remain stationary. the angle increased, greater ve- 
locity required hold the particle stationary move itup the plane. How- 
ever,as increases, the component the apparent weight normal the pipe 
wall that opposes the lift force decreases. Therefore some critical 
angle should possible have the lift force equal that any slight 
increase the value will cause the particle rise off the bottom the 
pipe. Since can readily measured, the lift force can determined for 
this condition impending motion. Alsofor this condition impending motion 
the drag force equal Thus both the lift and the drag 
force can determined. This method determining these forces has the ad- 
vantage that attachments need made the particle hold position. 

However, should noted that motion impends the direction the 
lift force both the normal force and the frictional force approach zero. 
Therefore order that rotation take place, indicates that the product 
approaches zero the lift force approaches thus questionable 
whether impending motion the direction the lift force will occur without 
rotation. However, actual observations indicated that particle can made 
rise off the bottom without rotation indicating that the necessary conditions 
were Satisfied, least for the tests run this investigation. 


EXPERIMENTAL PROCEDURES 


The experimental set-up used measure the lift and drag forces the 
procedure outlined was quite simple. Itconsisted primarily constant-head 
tank connected glass tube which could inclined various angles. The 
angles could measured within 1.0°. Parallel needle valves were used 
control the rate flow through the system. The fluid temperature remained 
relatively constant for all tests, and the fluid used all tests was water. 
schematic diagram the equipment shown Fig. 

was expected that for turbulent flow the pipe the motion the particle 
would erratic and measurements difficult obtain. For this reason was 
desired run all tests under laminar-flow conditions the pipe. Therefore 
particle low specific weight was needed that the desired conditions 
could obtained with relatively low velocities. meet this requirement 
nylon spheres with specific weight 1.1923 were used. The 
test facility was arranged that various sizes glass tubing could used. 
order that the flow the measurement station fully developed laminar 
flow, this station was located distance approximately from the en- 
trance and from the discharge end the tube. commonly used equation 
for estimating the transition length 


which the tube radius and the Reynolds number based the tube 
radius.4 The maximum transition length that occurred during the tests, 


“Applied Hydro- and Aeromechanics,” Prandtl and Tietjens, Dover Pub- 
lications, Inc., New York, 1957, 22. 
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estimated from Eq. was 3.3 ft. was therefore assumed that for all the 
tests the flow was fully developed the measurement station. Table sum- 
marizes the combination particle and tube sizes used the experimental 
investigation. 
For given test the angle fixed and the motion the particle the 
tube was controlled the rate flow from the constant-head tank. ex- 
pected, atsome critical angle the particle could made rise off the bottom 
the tube with essentially motion parallel the longitudinal axis the 
tube. The particle would lift off the bottom, then the lift force would apparent- 
decrease, and the particle would that there was actually “hop- 
ping” motion. 
The procedure followed each test was adjust the angle and rate flow 
until the particle had impending motion normal the pipe wall and record 
the data for this condition. The rate flow the tube was obtained col- 
lecting the fluid discharged from the tube and determining the time re- 
quired collect given volume. 


FLEXIBLE 
TUBING 
DISCHARGE 
TANK 
FLEXIBLE 
TUBING 


FIG, APPARATUS 


One disadvantage this method measuring the fluid forces acting ona 
particle that certain amount personal judgment involved establish- 
ing the condition impending motion the particle. order have some 
index the reliability the resulting data, each test with particular com- 
bination particle and pipe size was repeated least seven times and the 
standard deviation the lift and drag coefficients was evaluated. The stand- 
ard deviation from the mean values the lift and drag coefficients was less 
than the mean values for all tests. indicated Table the range 
values for the (d/D)-ratio was from 0.282 0.635. The Reynolds number 
range, based the pipe diameter and mean velocity, was 360 1,117. 
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CYLINDRICAL TUBES 


TABLE 1.—SUMMARY EXPERIMENTAL RESULTS 


ANALYSIS EXPERIMENTAL DATA 


was assumed that the variables influencing the lift force are the mean 


velocity fluid density fluid viscosity particle diameter and pipe di- 
ameter that 


Eq. can written terms three dimensionless parameters 


where isthe cross-sectional area the particle and the constant has been 
arbitrarily introduced put the equation more conventionalform. There- 
fore the expression for the lift force 


The Reynolds number based the pipe diameter. course Reynolds 
number based the particle diameter could have been used. However was 
felt that the velocity distribution preceding the particle, which function 
the pipe Reynolds number, would have considerable influence the lift and 
drag coefficients and thus was used. However the results can readily 


interpreted terms the Reynolds number based particle diameter 
since 


0.399 0.321 0.803 166 416 
0.493 0.559 1.135 380 1.27 
0.519 0.721 1.389 538 1,037 2.76 
0.562 0.638 1.135 402 715 3.39 1.73 
0.598 0.480 0.803 215 360 5.02 2.58 
0.635 0.721 1.135 411 647 2,43 
where the lift coefficient given 


where 


Since the lift and drag coefficients can evaluated from the experimental 
data, the problem reduces one determining the unknown functions and 
Table contains all the data and Fig. plot the experimentally de- 
termined values and versus for the variousdiameter ratios. Due 
the limited number particle and tube sizes available and the manner 
which the data were obtained was not possible hold either d/D con- 
stant for series tests. was therefore difficult correlate the data. 
However examination Fig. suggested that both the lift and drag coeffi- 
cients were inversely proportional was further assumed that the 
equations for the lift and drag coefficients could written 


and 


apparent that these assumed forms are correct then Rp!/3 
and versus d/D should correlate all the data, and these plots were 
used indicate the validity Eqs. and 13. was found that for the range 
variables covered the investigation the data appear correlate satis- 
factorily shown Figs. The empirical equations for the drag and 
lift coefficients which fit the data shown these figures are 


2.57 
and 
2.28 


should emphasized that although the relatively simple form for corre- 
lating thedata appears reasonably satisfactory, not likely that will 


; 
Similarly the drag force can expressed 
d 
| 
nail 


valid over the entire range Reynolds numbers that commonly inter- 
est. this investigation the range Reynolds numbers based particle di- 
ameter was 166 538. well knownfrom conventional drag-coefficient 
tests for spheres essentially unbounded fluid that this range Reynolds 
numbers represents transition region and that simple valid expression for 
the drag coefficient this range cannot extended into the regions where 
viscous forces are predominant into regions where viscous forces are neg- 
ligible. 


4.0 
0.500 
~ 
3.0 


0.635 


0.562 
~ 


FIG, 3.—VARIATION COEFFICIENTS WITH REYNOLDS NUMBER 


Figs. and indicate that for given Reynolds number both the lift and 
drag forces increase with increasing values d/D given value d/D 
both the lift and drag forces decrease the Reynolds number increases. 
These results for the drag coefficients are qualitative agreement with those 
obtained® McNown and Newlin their study the drag force 
acting spheres located the center tube. Data obtained 


“Drag Spheres within Cylindrical Boundaries,” McNown and New- 
lin, Proceedings, First Natl. Congress Applied Mechanics, 1951, 801. 
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from curves given McNown and Newlin are plotted Figs.6 and Itis 
noted that the drag coefficients for spheres located near the wall tube 
which the velocity approach has parabolic distribution are considerably 
higher than those for spheres located the center the tube which the 
expressions for the drag force acting sphere range placed 
the center cylindrical tube for both the case with uniform approach ve- 
locity and with approach velocity with parabolic distribution. They found 
that the drag coefficients for the case with the parabolic distribution were 
roughly twice those for the case with the uniform distribution the character- 
istic velocity mean velocity. indrag coefficients 
the same order magnitude those indicated Figs. and 
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From Figs. and apparent that the ratio the drag the lift force 
function only d/D. The plot this ratio versus d/D shown Fig. 


CONCLUSIONS 


this study method was developed determine the lift and drag forces 
acting spherical particle resting bottom circular tube. num- 
ber tests was performed confirm the feasibility using this method and 
determine the values the lift and drag forces. concluded that the 


“Motion Sphere and Fluid Cylindrical Tube,” Happel and Bryne, 
Industrial and Engineering Chemistry, Vol. 46, 1954, 1181. 
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proposed technique can used measure these forces accurately. addi- 
tion, Fig. shows that the lift force the same order magnitude the 
drag force and that the ratio drag although function thed/D-ratio, 

The range values parameters studied this investigation limited, 
but the same type equipment and technique can used extend the range 
with the proper choice fluids and particles. Probably the most important 
result this study the establishment the fact that for spheres near the 
wall tube alift force with magnitude comparable that of.the drag force 
does actually exist. thus apparent that the lift force should not over- 
looked studies related the incipient motion particles resting stream 
beds pipe walls. 


ACKNOWLEDGMENTS 


This investigation was conducted part the research program the 
Nuclear Engineering Group the Ames Laboratory. The assistance 
James Melcher conducting the tests and performing computations 
gratefully acknowledged. 


tye 
— 
a 
3 
|) 
ey 
SD 
; 
: 
oar 
A 


a 


June, 1960 


Journal the 


HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


MODELS PRIMARILY DEPENDENT THE REYNOLDS NUMBER 


SYNOPSIS 


The design, construction, and operation models closed-conduit fluid 
systems, and interpretation test data are discussed. Viscous and inertia 
forces the flowing fluids such models, and hence Reynolds number, are 
predominant factors. General rules for model size, construction, instrumen- 
tation, operating Reynolds number ranges, and data analysis are given. 


INTRODUCTION 


Hydraulic studies models closed-conduit flow systems provide re- 
liable means for obtaining useful and accurate design data. However, there 
are pitfalls the use models that must recognized and avoided the 
necessary accuracy and reliability are obtained. Thus, knowledge 
the similitude relationships between models and their full-sized counterparts 
essential. excellent discussion these similitude factors has been pre- 
account the details will not given here. will order, however, 
define what meant closed-conduit flow and state briefly the fundamental 
model relationships that must observed. 

Note.—Discussion open until November 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. June, 1960. 

Hydr. Engr., Div. Engrg. Labs., Bur. Reclamation, Denver, Colo. 


“Significance and Application Froude and Reynolds Numbers Criteria for 


Similitude,” Liu and Albertson, presneted the June 1959 ASCE Con- 
vention Fort Collins, Colo. 


“Engineering Hydraulics,” Rouse, John Wiley and Sons, Inc., New York, 1950, 
Chapter 
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Closed-conduit flow definedas flow enclosed system where free 
water surface exists, flow around deeply submerged body effect 
surface negligible. Typical examples structures involving 
closed-conduit flow include piping systems, pressure tunnels, control gates, 
pumps, and turbines. 

studies these systems, gravitational forces are not factor affecting 
the flow because direction always will down the slope the energy 
grade line and independent the slope the pipe itself. Surface wave motion 
and surface tension are nonexistent because there free surface, and, 
therefore, are not consideration. Elastic compression the fluid, which 
the case gases can significant, usually negligible with liquid flow. Thus, 
all the major forces and fluid characteristics affecting fluid flow, viscosity 
remains the predominant factor the usual closed conduit, hydraulic flow 
problem. 

The Reynolds number, the ratio the fluid inertia forces 
the fluid viscous forces. The number dimensionless, and the fluid 
density, denotes the flow velocity, characteristic dimension (usually 
the diameter pipes), and and are the absolute and kinematic fluid vis- 
cosities, respectively. Values for for fluids normally encountered tem- 
peratures are given Fig. 

Usually the Reynolds number prototype structure very large. Dupli- 
cating this number scale models impractical, particularly the same 
fluid usedon both structures, the size decreases, velocities must in- 
crease proportionately. With high initial prototype velocities, the problem 
rapidly would seem become impossible. 

Study the Reynolds number equation shows that the effects viscosity 
are particularly important when flow velocities and, hence, inertia forces, are 
low. But when flow velocities are high, inertial forces assume such great sig- 
nificance that viscosity becomes relatively minor. This basic factor that 
makes possible achieve reasonable similarity between model and its 
larger-sized prototype. making the model large enough, using ve- 
locities high enough, both, inertial forces can made predominate 
over viscous forces and reasonable similarity with the prototype will exist. 
The other requirement that the model large enough insure flow the 
fully turbulent range whenever full turbulence occurs the prototype. Thus, 
the procedure followed making usual model studies elements closed- 
conduit hydraulic systems build reasonably large, very accurate models, 
and operate them with flow velocities equal higher than scaled veloci- 
ties. Full turbulence generally encountered prototype structures, and the 
discussion this paper will limited these conditions. 


MODEL DESIGN 


first consideration designing model the selection the testing 
fluid. Water, the most generally available and perhaps most satisfactory li- 
quid, commonly used. Water also usually the fluid used prototype 
structures and naturally associated with model studies. many types closed- 
conduit testing, air can used the testing fluid with real savings time 
and costs.4 For this reason, low-velocity air tests should considered the 


Tests Using Low Velocity Air,” Ball, Transactions, ASCE, 1952. 
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FIG, 1.—KINEMATIC VISCOSITY FLUIDS VARIOUS TEMPERATURES 
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planning stages any model study. Sometimes may desirable take 
advantage fluids making rapid, low-cost air the early 
stages the study, and using water for the final tests and calibrations. 

The model should made large practicable, considering construc- 
tion and operating costs and the benefits that may derived. Increasing the 
size model usually enhances its usefulness and improves its accuracy. 
some point, however, the increase cost and the difficulty operation will 
offset any size advantage. Current practice follow precedent when avail- 
able and erron the side largeness insofar space and water supply 
mit. most closed-conduit models will fall within the scale ratio 
range 1:5 1:30. Undistorted scale modeling the rule. 

For reasonable similitude between model and prototype, usually 
sary operate the model Reynolds numbers approximately 106 


more. these values the frictional coefficient the Darcy-Weisback 
mula 


constant for manyof the most frequently encountered pipe roughnesses, and 
nearly constant for very smooth pipes (Fig. 2). the case valve conduit 
models using scaled heads and velocities, this requires about 6-in. pipe di- 
ameter. Consideration must given the whe- 
ther adequate discharge and head will available for the scale ratios being 
considered. this regard, must remembered that part the tests may 
made reduced gate openings and that heads much larger than scaled may 
necessary maintain reasonable values Reynolds number. 

The extent the total prototype structure represented the model 
must determined. Usually, interest centered around particular part 
system, rather than the system whole. such cases, necessary 
model only that particular portion the structure, with enough inlet and 
outlet piping produce flow conditions similar those the prototype. For 
example, the case control gate the midpoint long conduit dam 
approximately model conduit upstream from the gate, the gate 
itself, and about diameters conduit downstream might suffice (Fig. 3). 
the flow distribution were severely distorted multiple bends partially 
opened valve the start the conduit, longer upstream conduit would 
needed. another example, several control gates might connected 
branching manifold (Fig. 4). this instance, the characteristics the flow 
reaching each gate would affected the configuration the manifold and 
the number and location gates operating. Thus, would necessary 
model the entire manifold and all the gates represent fully conditions that 
could occur the prototype. studies were made high-head conduit in- 
let, might necessary provide pressure tank represent the 
voir (Fig. 5). studies with lower head box may adequate (Fig. 6). 
Each structure modeled must examined, and engineering judgment 
used, determine the model limits. 

Frequently, will advantageous model critical area only within 
structure and thison enlarged scale. example would gate slot 
for restraining the leaf control gate. relatively large slot can formed 
the wall test conduit, and accurate tests can made care taken 
have the opposite wall far enough away that there hydraulic 
ence (Fig. 7). Another example would the use pie-shaped sectors 
cular control valves (Fig. 8). Sectors less than 180° are the most common. 
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Time and unnecessary expense may saved selecting the scale ratio 
accommodate standard sized pipes and other model components that are al- 
ready hand. Odd scale ratios, such 1:17.35 willin way affect the accu- 
racy the model complicate the data analysis. 

Materials ordinarily usedin water models include bronze castings, plate and 
sheet brass, galvanized sheet steel, plastics, concrete, wax, and modeling clay. 
air models, wood and plaster are also satisfactory. When moderately high 
heads and velocities occur hydraulic models, rigid construction must 
used. Models circular cross section, such needle valves, use bronze cast- 
ings toadvantage. When the shapes are rectangular cross section, machined 


FIG, 9.—TRANSPARENT PLASTIC MODEL 


brass plates, doweled and bolted together, work well. advantage this 
build-up design that model changes can made readily. Transparent plastic 
used frequently when the flow within the structure must seen. Models 
entirely transparent are fairly common (Fig. 9). any case, the materials 
selected the design stages should combine durability, dimensional stability, 
ease construction, ease modification, and corrosion resistance. 

The inclusion measuring equipment such piezometers important 
design consideration. Often the requirement getting the instruments place 
within the model dictates the type construction used. Piezometers are rela- 
tively cheap when included the initial model construction, and should used 
generously offset oversights defining critical areas. Thought should 
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therefore given the necessary instrumentation and how will included 
the model the desired information will obtained. 

The testing fluid must supplied the model manner representative 
prototype conditions. Appropriate conduit length must provided between 
any supply system bends controls, and the model itself, that the flow dis- 
tribution will developed the proper extent the time reaches the test 
section. Straightening vanes are useful reducing the for spiralling 
flow within the pipes. 

The testing fluid may supplied directly from pumps, from surge 
constant head tanks. Ideally, would desirable have each model supplied 
through constant-head tank, but many cases this unnecessary and excel- 
lent results are obtained with the flow supplied directly from pumps. The 
characteristics the pumps being used have considerable bearing whether 
not constant-head tank necessary. the pumps are used directly, 
best run them appreciable heads and throttle the flow the model. 
this manner, better flow stability obtained. 

The waste from hydraulic models must either returned laboratory 
reservoirs the rivers lakes from which was drawn. This may 
done the most convenient manner and has little effect the model operation 
proper outlet conditions are maintainedat the model test section. Discharg- 
from air models usually present problems. 

Accurate and complete drawings are essential for expediting model con- 
struction and preventing time consuming, expensive errors. Drawings should 
contain sufficient informationto enable craftsmen build con- 
forms the designers’ specifications. Special drawings details are often 
necessary. All informationfor the installationof instruments suchas piezome- 
ters, pressure cells, velocity measuring devices, etc., should included, 
thereby obviating extensive changes after the model completed. fortunate 
cases where craftsmenare experienced hydraulic model construction, many 


structural details may omitted and much the construction method left 
the shop foreman. 


MODEL CONSTRUCTION 


Extreme essential for obtaining goodresults from 
any hydraulic air model. Care must taken that the materials, methods, 
and satisfactory and that tolerance checks are made frequently. 
Greatest accuracy should maintained where abrupt changes flow direction 
occur, and where flow velocities are high pressures are low. The stage 
the testing program affects the type and accuracy construction. early 
tests, during which many schemes may tried determine feasibility, con- 
struction need not careful the final stages, when data are used 
describe prototype operating characteristics. 

hydraulic models high head gates and valves, relatively thick metal 
plastic sections are usedand all important flow surfaces are machined. Thus, 
close tolerances and accurate component alinement can obtained. air 
models, wood, sheet metal, wax, and lighter plastic sections are permissible 
because the internal forces and the weight that the models sustain are small. 
The need for accuracy remains unchanged. For most water models made 
the Bureau Reclamation, United States Department the Interior (USBR), 
transitions from one shape another, manifolds, bends, etc., are 
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formed from gage sheet metal, 1/10 1/4 in. thick, hot-worked 
sheet plastic. cases where flow conditions within the system are ob- 
served, the transparent plastic sections are ideal. They are usually compara- 
tively expensive because wooden patterns must often made first, and the 
heated plastic formed over them, cooled, trimmed, joined, and flanged into 
final condition. The handwork increases costs, and the advantage visual in- 
spection must balanced against this extra expense. 

The inclusion piezometers often complicates the model construction and 
may restrict the choice materials. the case castings heavy metal 
sections, piezometer holes may drilled through the bodies and connecting 
taps can soldered intothe exterior surfaces. When number 
piezometers are grouped within limited area, may necessary touse de- 
vious passages prevent one piezometer lead from interfering with another. 
sheet metal construction, metal tubing inserted through the metal walls, 
soldered place with the tubing normal the wall, and then smoothed flush 
with the flow surface. All burrs must removed from the interior the 
piezometer opening.5 When concrete, plaster, wood, wax used for form- 
ing intricate flow surfaces, desirable embedin the material metal tubes 
that extend the surface, that the piezometers will have durable, smooth, 
square-edged openings. When materials are suitable and superior craftsman- 
shipis available that damaging donot occur, holes can drilled 
directly into the material. Usually the should larger 
inside diameter and should extend least diameters from the opening 
without change incross section alinement. Piezometers smaller than 1/16- 
in. diameter may used, and may desirable curving surfaces where 
many are crowded into small area, but they are prone becoming plugged 
with small bits material, and tend induce damping water models that 
reduces the sensitivity measurement. The smaller piezometers work well 
air models. 

Flexible tubing about 1/4-in. inside diameter usually used connect the 
piezometer taps the pressure measuring instruments, usually water man- 
ometers. Rubber tubing has been used extensively, but relatively short-lived 
and has anopaque quality that prohibits observing interior conditions. Recent- 
ly, flexible, transparent plastic tubings have become available and are highly 
recommended. advantage the transparent tubing that entrapped air 
bubbles can detected easily and removed from the line. 

When measurements are made rapidly fluctuating pressures, care 
must matching piezometer diameter, pressure lead diameter, pres- 
sure lead length, the pressure cell itself, and the Failure 
this may result obtaining records that reflect serious distortions the ac- 
tual pressure history. Deductions made all innocence and confidence from 
such records would, course, erroneous and misleading. general, 
should remembered that the pressure leads should short possible, 
rigid, and least inside diameter. All and gas 
bubbles should bled from the system, and the temperature the cell should 
not vary appreciably during the calibration and tests. 

“Piezometer Investigation,” Allen and Hooper, Transactions, ASME, 
Hyd 54-1. 

Note the Evaluation Design Transducer for the Measurement 


Dynamic Pressures Liquid Systems,” Barton, Instrument Notes, Statham 
Labs., No. 27, June, 


“Attenuation Oscillatory Pressures Instrument Lines,” Iberall, Re- 
search Paper 2115, Natl. Bur. Standards, Vol. 45, July, 1950. 
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Provisions should made for installing any other measuring equipment 
needed the model. For instance, may desirable make velocity tra- 
verses across the conduits, and pitot tubes hot wire anemometers may have 
inserted through the walls. other cases, the injection dyes fine 
streams air bubbles may used show flow patterns. 

Flow-measuring devices such Venturi meters, orifice meters, weir 
boxes are required for determining the rate flow. Much depends the ac- 
curacy these flow measuring devices and carefully calibrated, properly placed 
equipment essential. Quality metering devices are available and the se- 
lections can based accuracy desired, cost, convenience use, and ex- 
pected service life. 

Properly fitted, tightly connected pipes, fittings, and transitions should 
used between the the flow-measuring device insure that all mea- 
sured flow reaches the model test section. This requirement applies equally 
water and air models. Leaks air models are often difficult detect be- 
cause they leave visible trace. water models, outside leakage and drips 
are easily detected, but leakage air into zones subatmospheric pressure 
less obvious. 

assembling the model components accurate alinement must maintained. 
Adequate supports, provided the model, should placed provide 
minimum interference when alterations the model are made. Usually 


convenient place the model waist height for ease construction, oper- 
ation, and observation. 


MODEL OPERATION 


program must carefully planned make full use the model. 
are necessary reveal defects, leaks, other inadequacies. 
First tests whether not minor changes are needed make the 
model and the instrumentation better fitted for obtaining the desired data. This 
phase should not hurried because economically important certain 
the model really represents prototype conditions and that the instrumentation 
satisfactory. during this adjustment phase that the experimenter be- 
comes acquainted with the peculiarities the model andadept its operation. 
Testing should consist examining systematically each proposed, feasible 
design for possibilities for improvement flow conditions, cost reduction, and 
ease maintenance. The experimenter must exercise patience, imagination, 
and ingenuity and capable interpreting the model results correctly. The 
data should analyzed concurrently with the testing prevent accumulation 
unnecessary data, and permit intelligent decisions concerning succeeding 
steps the testing program. relatively complete and detailed diary the 
study important. This diary should include dates, general notes 
mance, changes made, the model data, and photographic records. Often ap- 
preciable period time will elapse between the time model study and the 
time when report prepared. Without record such the diary provides, 
many factors could forgotten overlooked. 

Measurements pressures, discharges, velocities, etc., must made 
with precision. Much depends the values recorded the data sheets, and 
upon pertinent remarks such pressures being particularly unsteady cer- 
tain locations. The type pressure and velocity records made depends upon 
the use intended for the data. many cases, small pressure fluctuations are 
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not significant and the average pressures are the ones desired. Damping de- 
vices may used these pressure lines make the reading simpler. Other 
cases willoccur where transient pressures are desired. These measurements 
must made with special equipment such pressure cells and electronic 
recorders that are capable following the pressure swings. appreciable 
damping should occur during these measurements. all pressure readings, 
datum points gage zeroes must recorded that actual pressure values 
the piezometers will known. The testing must extended over the en- 
tire range equivalent discharges and operating conditions that will occur 
the prototype. Frequently desirable extend the range the tests 
heads much higher than the equivalent heads anticipated for the prototype 
fully exploit the model. 

addition routine testing immediate interest, the experimenter should 
alert obtain general information that will value later date. Re- 
sults experimentson many models are required for assembling general de- 
sign information. Time and money may wasted all the useful information 
not obtained from each model under test. 

Studies made using air flowing fluid usually require that the flow ve- 
locities maintained below approximately 300 fps. higher velocities are 
used, compressibility would have considered. The density the air will 
vary from day day with changes intemperature and/or barometric pressure. 
Thus, pressures are read water manometers feet, and are converted 
feet flowing fluid (air), the density air that water must 
known. The air temperature and barometric pressure must therefore re- 
corded during air studies. 


INTERPRETATION RESULTS 


care taken operate the model Reynolds numbers about 106 
more, the data may analyzedusing the same relations for models based 
the Froude law. Thus, direct geometric linear values possible. 
the model were run with scaled heads and discharges, the model pressures 
would multiplied the scale ratio obtain prototype pressures. Time 
intervals and flow velocities are obtained multiplying model values the 
square root the scale ratio. The discharge obtained multiplying the 
model flow the five-halves power the scale ratio. 

the modelis tested heads and discharges that differ from scaled proto- 
type values, the data may placed dimensionless form and then applied 
prototype conditions. example shown Fig. 10. typical dimension- 
less pressure factor relation 


where the head the piezometer question, the head refer- 
ence piezometer, and the Q/A velocity the reference station. The 
reference station usually taken convenient point just upstream from the 
test section itself. Other pressure factor relations may used conditions 
dictate, provided they involve pressure relationship from base point 
parameter, such the velocity head, that describes the flow. 
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Data the foregoing dimensionless form applied follows: Determine 
from the relation 


the rate flow the prototype structure. Eq. significant area 
the gate conduit, denotes the coefficient discharge, and the 
total head differential across the structure. After obtaining compute the 
velocity head the reference station. The pressure factor value for the par- 
ticular piezometer then multiplied this velocity head give the pressure 
change feet water piezometer the reference station. The pres- 
sure the reference station found adding any back the struc- 
ture the head caused the presence the structure, and given the re- 


lation The computed pressure change, added algebraically 


the reference pressure, gives the prototype piezometric pressure expect- 
ed. Data obtained from air model tests are best expressed pressure factor 
form, using care that all values are expressed feet the flowing fluid. The 
pressure factors are dimensionless and are applicable directly prototype 
conditions using other fluids. 

Often negative, subatmospheric, pressure heads are obtained from models 
which, when scaled directly, yield values too low possible the proto- 
type. example would 3-ft subatmospheric pressure head measured 
1:20 scale model. Direct scaling yields prototype head below at- 
mospheric pressure. This unrealistic because pressures lower than ab- 
solute zero can obtained. Atsea level, this about -34 ftof water. Actually, 
water begins when the pressure reduced vapor pressure, about 0.7 
absolute, and cavitation occurs. Thus, any scaled-up model pressures that 
are subatmospheric the extent about indicate that cavitation 
willoccur the prototype. cavitation indicated over appreci- 
able areas the structure, large-scale separation likely and the flow pat- 
tern may materially changed. Specialized testing necessary evaluate 
such structures. the usual case, the mere suggestion pressures near the 
cavitation range sufficient justify redesign that noise, vibration, and 
damage caused cavitation will not occur. 

some cases, negative pressures limited extent may acceptable 
prototype structures. When flow velocities and pressure fluctuations are moder 
ate small, indicated pressures below atmospheric are not considered 
objectionable. usually dangerous allow pressures lower than this be- 
cause unexpected surface roughnesses, vorticity, and flow turbulences may 
momentarily lower local pressures the cavitation range. 

Studies cavitation potentialon particular structures are frequently made 
and often require special test equipment and instrumentation because the 
complexity the problems. Certain parameters are commonly used express 
cavitation potential structure. One the cavitation index 
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which the pressure reference point, the vapor pressure 
relative the atmosphere (negative), and the velocity reference 


point. Another index especially useful for systems such valves pipelines 


tance D), hyp the vapor pressure, and the total head upstream. 
This parameter has the advantage describing pressure conditions which 
must exist with given flow velocities and valve characteristics cavitation 
With either parameter, for the same value model and 
prototype, the pattern cavitation will the same. But since frequency 
vibration and rate pitting vary with velocity, care must used predict- 
ing prototype structural behavior.9 

The head loss through the system obtained taking pressure and velocity 
measurements appropriate points the system. Similarly, the coefficient 
dischargefor the structure may obtained with appropriate measurements 


pressure, discharge, and valve gate opening. This coefficient common- 
expressed 


which discharge, head drop across the structure, and de- 
fining area similar that the conduit, the control body, the opening be- 
neath the gate leaf. Force measurements, such downpull gate leaf, 
torque turning vanes, thrust walls, may obtained applying mea- 
sured pressure heads areas that will affected, making direct force 
measurements. the latter method used, care must taken eliminate 
unwanted model friction and extraneous vibration. 

The usefulness air models studying closed conduit systems has been 
discussed. Such models discharge into the atmosphere and therefore operate 
“submerged,” which analogous water model discharging under water into 
large pool. Air-model data must therefore interpretedas submerged flow. 
This closed conduit systems, but must recognized and ac- 
counted for model studies are made prototype structure such gate 
that discharges water freely the 

Hydraulic models sometimes combine both open and closed conduit flow 
conditions. outlet works illustrates case where studies are made the 
control valve gate and the stilling basin downstream. Both the open and 
closed conduit portions the model are analyzed the Froude relationships, 
the operating Reynolds numbers are the appropriate range. 

difficult build model with equivalent dynamic response that 
the prototype. This does not preclude making studies the stability the 
system reasonable predictions whether not resonance and vibration 


“Cavitation Characteristics Gate Valves and Globe Valves Used Flow Reg- 
ulations under Heads about 125 Feet,” Ball, Transactions, ASME, Au- 
gust, 1957. 

“Cavitation Scale Effect,” Knapp, Proceedings, Seventh Genl. Meeting, 
A.H.R., Lisbon, 1947, A6. 

“Air Model Studies Hydraulic Downpull Large Gates,” Simmons, 
Jr., Proceedings, ASCE, Vol. 85, No. January, 1959. 
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will occur. Simultaneous recordings measurements appro- 
priate points within system can show pulsations and/or surge will pres- 
ent. These data may applied mathematically the prototype structure for 
analysis possible resonance. such tests, particular care must taken 
supply water the model without pump other mechanically induced 
sure fluctuations. records will requiredon the incoming pipe- 
line show whether not such pressure fluctuations exist. Special and de- 
tailed data analysis are required such fluctuations are entering the system 
complicate the pressure conditions. 

Graphs, charts, and other plots are usually needed analyze and illustrate 
the data obtained. Final results find more usefulness when expressed di- 
mensionless parameters that may applied other geometrically similar 
structures. Plots showing specific pressure conditions and other 
istics applicable the structure under study are also made. 

Reports constitute the chief instrument through which the benefits the 
model study are transmitted designing engineers, and the profession 
general. They must complete, accurate, timely, and the point. However, 


brevity virtue attained the expense completeness, clearness, 
accuracy. 


CONCLUSIONS 


understanding Reynolds similitude factors, and proper application 
them designing, operating, and interpreting closed conduit hydraulic models 
will lead reliable and accurate predictions prototype performance. 

Models used the usual studies hydraulic control elements and con- 
duit systems need not excessively large, operated with extreme flow ve- 
locities yield reasonable model-prototype similitude. rule thumb, 

model operating near Reynolds number about 106 adequate. 
Surfaces model flow passages must usually smoother than the 
prototype surfaces for friction similar. short portions 
conduit systems, full similitude this respect not important, and good re- 
sults will obtained the velocity distribution entering the section turbu- 
lent and fully developed. models with long water passages, more effort may 
required obtain friction similitude, required, because friction would 
relatively great importance. 
Particular care must taken produce accurately shaped flow sur- 
faces, and maintain precise alinement the model components. 

The selection, location, construction and maintenance measuring de- 
vices, such piezometers and recording instruments, must done 
and knowledge the type data obtained. 

Tests should continued the point obtaining all feasible information 


from each model that the time and expense the model construction will 
exploited fully. 
Detailed diaries logs should kept the studies, and comprehensive 
reports prepared make the material available the engineering profession. 
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HYDRAULIC MODELS THE ST. LAWRENCE POWER 


BRYCE.!—The author wishes express his appreciation the dis- 
Brown, both which contained much additional useful information. men- 
tioned previously, the comprehensive nature the paper necessitated gen- 
eral coverage the subject and space did not permit detail the 
various problems encountered. However, both discussions, interest was 
shown the designof energy dissipation works and more information was re- 
quested connection with the investigation hydraulic loading gates. 

With regard the energy dissipation the Iroquois Control Dam, Messrs. 
Fenwick and Brown only the sluices were thus protected. 
they surmise, very substantial economies were effected this treat- 
ment. careful analysis supplemented model tests indicated that only dur- 
ing construction would energy dissipator required, and this high head 
period the flow could confined sluices. This was later corroborated 
actual operation. These discussers also note that the pattern gate openings 
can have considerable effect energy dissipation and the magnitude 
the currents downstream. This was forcibly illustrated Iroquois during con- 
structionand even after the power pool was raised. Initially during these peri- 
ods, fairly severe currents were observed below the dam which proved trouble- 
some navigation. intensive series model tests was run the 
Iroquois Dam model, and patternof gate openings was developed maximize 
energy dissipation and minimise downstream currents. When these patterns 
were put operationat the dam remarkable improvement occurred and con- 
ditions acceptable navigation were produced. 

Messrs. Fenwick and Brown raised the question thenecessity 40- 
wide end sill for the Long Sault Dam dissipator and whether longer apron 
and shorter sill might improve conditions. mentioned the paper, over 
thirty different designs were tested which included several combinations 
the above variables, and was considered that the design selected was the best 
view all the circumstances. Some the factors which influenced the se- 
lection were: the depth the bucket was limited from geological 
ations; was desired have the apron shortas possible, consistent with 
satisfactory performance; and the issuing velocities were still rather 
high, the end sill was made wider than hydraulic conditions alone would require 
effort ensure that, undercutting did take place, sufficient sill would 
remain function properly. 

Both Mr. Berry and Messrs. Fenwick and Brown requested further details 
regarding the investigation hydraulic loading gates. Although space will 

Hydr. Engr., Hydro-Elec. Power Ontario, Toronto, Canada. 
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June, 1960 


not permit complete discussion this subject, the writer feels that certain 
information additional that given the paper might prove useful. Accord- 
ingly, and with reference the original paper, the following remarks are ap- 
pended. 

Port Gates Long Sault Dam.—The results the tests these gates are 
given Fig. the paper. The gates were ft, in. high, wide, and 
in. deep, and had the skinplate the downstream side. The tunnel ports 
were high the upstream and highat the downstream end and were 
wide. Each gate was capable shutting off maximum flow through 
the port 14,000 cfs travelling average velocity fps and was 
fully submerged operation. All tests the gates out head 
the sill and with free discharge conditions the outlet. 

Two points not mentioned the paper might prove interest. The first 
that when tests were carried out the original gate design, the flat bottom 
caused the spring point the flow move from the downstream edge towards 
the upstream side under certain conditions. When this took place oscillat- 
ing vibrationof impressive proportions occurred. This behaviour agree- 
ment with that cited Messrs. Fenwick and Brown their discussion. The 
the effect the downpullof recess made the upstream face 
the dam above eachport and immediately each gate when inthe wide 
open position. The recess eachcase was wide, high and deep, 
and was found reduce the downpull tons. believed that the re- 
duction downpull the recess due balancing the pressure across the 
top seal. However, the recess was extended beyond the top the gate the 
fully open position, there was tendency for the gate move outwards from 
the dam. The significance the clearance the downstream side the gate 
was also mentioned Messrs. Fenwick and Brown. 

Following the tests reported the paper, additional investigation was 
made determine the effect several variables. brief was found that, 
without air venting, downpull varies directly with head, and that air venting 
increasingly effective reducing downpull the head increases. was also 
found that approach conditions have considerable effect the vertical forces 
the gate. For example, lowering the bed upstream some reduced the 
downpull over 20%. 

Powerhouse Headgates.— The headgates tested were those for the Canadian 
powerhouse which are somewhat different design than those for the United 
States powerhouse. Fig. shows the gate dimensions and the details one 
the three water passages leading the scroll case. The gate skinplate 
the downstream side and the gates operate fully submerged condition. 
was determined that the critical conditions for hydraulic loading the gate 
would occur under runaway conditions with two the inlets completely closed 
before the third gate started downward. The tests were performed with con- 
stant headwater and tailwater levels. The wicket gates pass 14,000 
cfs the beginning the test (runaway discharge) and locked position. 

Two the passages were closed leaving the gate under test toclose the re- 
maining opening which was found tobe carrying about 13,000 cfs fps when 
the test commenced. Fig. the results show that the initial downpull 
tons increased maximum 119 tons gate opening 8ft. This may 
compared toa gate weight tons. was considered that the gate bottom 


was sloped muchas structural design would permit and reduction load 
was found possible. 
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DISCUSSION 


Massena Intake Gates.—At the Massena Intake, two gates were tested, the 
emergency gate which had skinplate the downstream side, and the service 
gate which had the skinplate the upstream side. Fig. shows the arrange- 
ment and dimensions the gates and the details the water passages. The 
tests were carriedout constant headwater level and tailwater 
levels. The maximum test flow was 21,000 cfs producing velocities over 
fps through the port. 


WEIGHT TONS 


MAIN AIR VENT 


SECTION 


GATE 


HYDRAULIC TONS 


18.—MEASURED VERTICAL FORCES CANADIAN HEADGATES 
FOR THE ST. LAWRENCE POWER PROJECT 


Emergency Gate Tests.—The emergency gate upstream from the service 
gate and completely submerged operation. shown Fig. 19, maxi- 
mum downpull forces tons and tons were measured for the two heads 
tested. may noted that the fully open position downpull about 
tons was measured, about 140% the gate weight. Comparing these downpull 
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forces with those measured the Long Sault Dam port gates would appear 
that they are somewhat greater than would expected. believed that the 
two horizontal sections the bottom, which accommodated the wheel, axles, 
increased the downpull force. Also, this gate 50% higher than the Long Sault 
gate and the drag forces are probably greater. 

Service Gate the service gate similar size tothe emer- 
gency gate, and the test conditions were identical, the results the tests were 
greatly shown Fig. 19. Downpull forces and 


SECTIONS 


EMERGENCY 45° SLOPING BOTTOM 
GATE 
335 HEAD 
HEAD 
a 
oO 
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TONS 


FIG, HYDRAULIC FORCES MASSENA INTAKE GATES 


tons respectively with the two heads, but uplift forces were encountered 
and tons respectively. violent vibration was observed and near 
gate opening the higher tailwater level tested, when variations 
force plus minus tons were observed. This tendency was reduced but 
not eliminated the lower tailwater level. believed that the uplift was 
due the considerable back-pressure the gate from the tailwater, and that 
the vibration was due unstable conditions certain gate openings when 
rather violent pumping action took place the tailwater. 
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DISCUSSION 


Sluice Gates, Long Sault Crest and Diversion these 
gates are similar that wide, thickness and have bottoms 
which sloped upward from the skinplate the direction flow angle 
approximately 45°. All the gates had upstream skinplates and were unsubmerged 
operation. The Iroquois gates were sheathed also the downstream side 
permit heating for winter operation. Fig. the results all the tests 
and the gate dimensions and sketches the gate settings are given. should 
noted that the head referred that the sill for the Long Sault Crest 


VARIES 
FROM 
30' 
LONG SAULT DAM LONG SAULT DAM 
CREST GATE SLUICE GATE DIVERSION SLUICE GATE 
100 TONS WEIGHT TONS WEIGHT TONS 


HEAD 


NOT 


OPERATING HEAD 25FT 
220 FOR ALL GATES 
LOADS ARE DUE 
HYDRAULIC ACTION ONLY 
WITH BUOYANCY 
EXCLUDED. 

HYDRAULIC DOWNPULL 

TONS 


FIG, 20.—MEASURED VERTICAL FORCES CREST AND SLUICE 
GATES FOR LONG SAULT DAM AND IROQUOIS DAM. 


gate, and the differential head theother two cases. was concluded from 
the results that the case the Long Sault Crest gates and the Iroquois gates, 
the downpull was caused drag forces only, atmospheric pressure oc- 
both gates the point maximum downpull. The increased 
downpull the Long Sault Diversion gates was believed due the restricted 
setting great mass water was thrown into the exposed 
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side the gate, thereby increasing the load. This was confirmed testing 
this gate the Iroquois setting where the downpull was greatly reduced. 

obtain some information the effects bottom design gates this 
type the Iroquois gate was tested with flat bottom. the tests, the point 
where submergence the gate bottom began, violent vibration occurred. 
The maximum downpull was also increased from tons tons. 

While the tests previously described means cover the subject ver- 
tical forces gates, would appear that certain general conclusions may 
drawn. Vertical hydraulic forces have been shown very large num- 
ber gate installations, particularly when the gate has downstream skin- 
plate and submerged. 

With downstream skinplate, downpull the dominant force which may 
many times exceed the gate weight. Substantial downpull forces may occur 
with the gate fully open suspended close submerged port. Uplift 
appears relatively minor with this type gate. The magnitude the 
downpull appears vary with the head the gate and the size the gate and 
greatly affected the shape the gate bottom, being greatest with flat 
bottom. While the major cause downpull undoubtedly the reduction 
pressure under the gate due the velocity, has been shown that drag forces 
and approach conditions also have appreciable effect. Effective air venting 
and equalizing pressures across projecting top seal appear reduce down- 
pull, but sloping the bottom members the gate downward the direction 
the flow produces the greatest reduction. 

With upstream skinplates, downpull not normally large, but uplift may 
important the gate submerged. Sloping the gate bottom upwards the 
direction flow appears beneficial reducing unstable operating conditions 
part gate openings. 
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GROUND-WATER PROBLEMS NEW YORK AND NEW 


Closure Joseph Upson 


JOSEPH UPSON.!—The discussions the writer’s paper Messrs. 
Biemondand Robert Thomas reveal twoquite different categories ground- 
water problems. Mr. Thomas’ comments elaborate the need for adequate 
investigation the geologic, hydrologic, and economic features pertinent 
the development ground-water resources. Solutions these general prob- 
lems are certainly fundamental the early stages planning water-supply 
installation, whether for irrigation, industry, domestic supply. 

contrast, Mr. Biemond points out much more specific type problem 
which lies the narrower realm the minute relationships between ground 
water, with its dissolved, suspended, and bacterial loads, and the small-scale 
environment its enclosing deposits. The effects these relationships some- 
times appear operation installation and can almost critical 
the amount available water. Both categories are applicable and pertinent 
development the ground-water resources New York and New England. 
Few workers with, and students of, water are knowledgeable the fields 
both categories. 

For the readers who may not fully appreciate Mr. Biemond’s final reference 
advantages artificial replenishment, the Amsterdam, The Netherlands, 
water works has been pumping water from the Rhine within The Netherlands 
and, with minor pretreatment, infiltrating into segment the ground-water 
reservoir beneath The Netherlands coastal dunes since 1957. Apparently there 
has thus far been clogging other difficulties. 

The writer extends thanks and appreciationto Messrs. Biemondand Thomas 
for interesting and informative discussions. 


June, 1959, Joseph Upson. 
Research Geologist, Ground Water Branch, Water Resources Div., U.S. Geol. Sur- 
vey, Mineola, 
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ENGINEERING APPRAISAL HYDROLOGIC 


Task Group Hydrologic Data the Committee Hydrology. 


DONALD VAN TUYL,! thorough explanation Mr. Kaz- 
mann yields ground water aquifers should greatly 
assist those considering the types hydrologic data reviewed the report 
the Task Group Hydrologic Data. Classifying “estimated yield” interpre- 
tive data was intended the distinction between quantities determined 
compiling other values and quantities determined more less routine 
measurement analysis. Primary emphasis was given throughout the report 
basic and analyzed data with only bare mention the interpretive classifi- 
cation. 

The twoterms which Mr. Kazmann discusses are usedin ground water com- 
putations without exact definition, clearly indicates. Both may found 
have specific values during engineering analyses, however, which may also 
involve economic factors. 

The term “interpretive data” the Task Group report should viewed 
intermediate product raw materials (data) but not endin itself—a value 
knowledge usedwith other items obtaining further answer. 
item datathat results from interpretation and extrapolation basic data 
and employed subsequent hydrologic applications, therefore, may 
classed interpretive under this concept. 

For example, computed flood stage for particular frequency the re- 
sult numerous compilations basic and analyzed data, but meaningful 
only when applied item data engineering problem. member 
the TaskGroup suggests that storage-yield curves for surface water develop- 
ment are also interpretive data. 

One reference inthe another publication, unfortunately, 
was out date when the report reached printed form. Reference No. should 
read: 


Langbein, and Hoyt, G., 1959, “Water Facts for the Nation’s 
Future,” The Ronald Press Co., New York; sponsored the Conser- 
vation Foundation. 


July, 1959, Donald Van Tuyl. 
Task Group Hydrologic Data, Committee Hydrology, and Water Re- 
sources Asst., Chamber Commerce S., Washington, 
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HYDRAULIC CHARACTERISTICS GATE SLOTS# 


writer found Mr. Ball’s paper most interesting, es- 
pecially concerning the large number model test results that has been 
able furnish various different type slots. 

certain number the results with the different type slots obtained 
Mr. Ball agree with those found during the detailed program investigations 
carried out during the last the laboratories the Societe Grenobloise 
d’Applications Hydrauliques, Grenoble. The investigations were 
especially concerned with the problem partial openings vertical lift gates 
under very high heads and the resulting cavitation conditions the gate slots. 

The writer would, however, like stress that absolutely essential that, 
order interpret thetest results obtained correctly, the similitude between 
the prototype and the model must established perfectly. 

this not done, number cases the interpretations may prove 
rather faulty. Very great care this account has taken designing 
the test installation, that the Froude similitude and the Reynolds number are 
sufficiently high. The test installation used the laboratory had the clear 
gate dimensions: 600 high and 260 wide. The gate model was hy- 
draulically operated servomotor control. The cavitation circuit was de- 
signed pass 1200 liters water per sec. With this type circuit having 
satisfactory pressure control, was possible obtain reasonably accurately 
the approach cavitation breakdown. Also, order that the results obtained 
are not faulty, essential that the circuit clear parasite cavitation 
bends, the pump, etc., that the observed changes the regime are the 
gate operation and hydraulic conditions only. Further, water with some air- 
content used, the circuit should provided reabsorber, that 
bubbles not collect dead pockets. 

The writer completely agrees with Mr. Ball that the slot problemis essen- 
tially detailed study pressure intensity and distribution. The heads, the 
dimensions well the downstream hydraulic conditions will determine 
whether there cavitation not. The intensity cavitation will the same 
time depend the gate opening. For different hydraulic and slot conditions, 
the critical openings vary. For number existing installations, the gates 
have not shown signs cavitation, because they have never functioned under 
the critical openings for the hydraulic conditions imposed. 

The writer feels that the correct mode attacking the problem was not 
find typical design which would fool-proof under all different operating and 
hydraulic conditions but study each particular case separately. 


Engr., Sogreah, Grenoble, France. 
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The series tests that have been carried out the writer different 
types slots and hydraulic conditions imposed have proved that for 
effective heads over with downstream undrowned (free discharge) and 
with typical classical vertical lift gates, freedom from cavitation cannot had 
under all possible openings. This condition also applicable for arrange- 
ments having upstream contraction and air duct tubes feeding along the height 
the slots. Actually, the best arrangement with upstream sealing and nor- 
mal slots without offset. With this disposition, would possible operate 
the gates without large cavitation risks for higher heads, but only for inter- 
mediate (between 1/3 and 2/3) and full gate openings. 

however the downstream jet drowned (backpressure), the results are 
most unfavorable whatever the position the sealing, the type the slot, 
the gate opening. 

the vertical lift gates are required function partially open under very 
high heads, then solution other than the classical gate has envisaged. 
During our investigations, was specified that solution should foundwhich 
should not only suitable from the hydraulic point view, but had feas- 
ible mechanically and from the economic stand point, the cost had not ex- 
ceed those the classical vertical lift gate. The solution baffles exten- 
sions bottom corner gate leaves (in any case not very healthy solution) 
very large offsets the slots had rejected. After considerable 
amount study and tests, gate design was evolved which satisfied the differ 
ent conditions imposed. 

brief mention this solution was made Messrs. Zaky and Sauvage 
Saint Marc the last Congress I.A.H.R. 

impossible give complete idea this solutionin this short note 
comments. Briefly speaking, the gate leaf knife edge designed up- 
stream the gate slot, being supported sort parabolic bridge struc- 
ture. The slots from the hydraulic point view are nearly uniform flow 
without large variations pressure along the height the slot. The results 
that were obtained during the tests were most encouraging. 

case partial gate openings, just downstream the gate, one has over 
the resulting jet surface either air (free surface flow) mass water 
(drowned jet) which the pressure distribution practically hydrostatic. 
case free surface downstream, the downstream (back pressure) the 
height the jet itself above the gate sill. case drowned jet, the down- 
stream the height water the air duct, also based the sill level the 
gate. 

obtain accurate results, very delicate pressure measurements have 
taken. exact idea pressure fluctuations essential. know the 
average pressures not enough. 

The cavitation coefficient can described as: 


The numerator the difference between the absolute pressure downstream 

the jet, which 

itself the difference between the pressure just the upstream and down- 

stream side the gate. The sill taken datum for all heads. 


the gate and the vapour pressure hy. The denominator 
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DISCUSSION 


The classical vertical lift gates with upstream sealing and normal slots 
without offset and downstream undrowned cavitate around 0.25 difficult 
openings (less than 0.2 and greater than 0.7) and 0.1 favorable open- 
ings. With downstream jet drowned, the rarely under 

The writer trusts that the comments give least rough idea the way 
the problem gate slots has been studied and that they will help compare 
with the results obtained Mr. Ball. 
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DISCHARGE FORMULA FOR STRAIGHT ALLUVIAL 


Mostafa, and James Culbertson and Carl Nordin, Jr. 


extremely interesting investigation, from the 
practical point view, because provides engineers with easy homogene- 
ous calculation method that covers nearly all the experimental results known 
today, which not the case with certain number earlier investigations. 
Like earlier publications, this article shows that the problem not only very 
complex, but also that despite the many experimental results that are known 
and applied, their number and organisation are still grossly insufficient. 

spite this well known weakness, that affects all, the authors—like 
many their predecessors—have given the temptation generalise the 
solution and give vaguely theoretical appearance. They have, for in- 
stance, considered parameters for which little experimental data, none 
shape, grain and fluid specific gravities, viscosity, gravity and temperature, 
none which have really been considered the tests under analysis. Some 
doubts therefore remain the reliability the results, these are pre- 
cisely the parameters applied. 

For instance, although the introduction the drag coefficient deduced 
from the fall velocity grain still water probably valid for the more 
less spherical grains used the tests, doubtful whether would still 
apply, if, say, pins were used instead. Under these conditions, the form effect 
isolated pin under free fall conditions would probably bear relationship 
that pin lying flat the surface the water. 


Considering sphere, also doubtful whether the sudden change the 
drag coefficient occurring when 


under free fall conditions incalm water also valid for bed spheres. The 
difference state the flows and the geometry the system too 
great. the writer’s view, the introduction the fall velocity isolated 
grain not representative theform effect the when forming 

further doubtful result, which probably due the fact that the density 
the grain has not been considered over wide range, that the apparent 
weight the grain water does not affect the velocity under turbulent con- 
ditions (Eq. 26). other words, for the same rate flow, the same slope, the 


November, 1959, Hsin-Kuan Liu and Shoi-Yean Hwang. 
Engr., Sogréah, Grenoble, France. 
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same canal, and grains the same diameter but widely varying specific 
gravities, (with high sediment load one case and none the other) the ve- 
locity supposed the same. this really so? could that this 
result was obtained because the tests considered were not carried out over suf- 
ficiently wide specific gravity and sediment load ranges 

the title the article suggests, considered that the results the 
investigation are very true for natural alluvial channels, which such parame- 
ters gravity, grain and liquid specific gravities, grain shape, viscosity and 
temperature only vary slightly. However, contrary the apparent possibili- 
ties shown the presentation the results, some doubts arise soonas one 
leaves the scope the basic experiments. According the author’s conclu- 
sions, should possible settle the question once and for all carrying 
out tests very wide scale. 

The writer feels, however, that such grain shape, specific gravity, 
sediment load and viscosity should investigated very thoroughly. Tests 
should with variety shapes like pins, discs, lenticular forms, 
and forth, with relative densities varying between 100, and with very 
high sediment loads. 

“La Houille Blanche,” (June, 1959). Rottner collected 2,500 experi- 


mental results and succeeded classifying them terms the following 
parameters: 


which the sediment load, and the remaining the arti- 
cle under discussion. 

Rottner’s investigations not only enables the same determination 
made, but also gives the sediment load. thus more complete and opens 


FIG. 1.—DATA FOR STRAIGHT ALLUVIAL CHANNELS 


DISCUSSION 


new prospects certain number peculiarities and zones far inade- 
quately investigated empirical methods. The graph obtained, Fig. shows 
certain peculiarities that will readily noted specialists the subject. 


mula for alluvial channels ingenious, but his basic assumption that various 
empirical can synthesized into one rational form does not ap- 
pear justified. This one universe which order must brought not 
through synthesis but through catharsis. 

The author used the Blasius formula for grain sizes smaller than 0.1 mm. 
The Blasius formula valid only the range Reynolds numbers less than 
80,000, which means that for 1.12 10-5 for 65° (as used Liu) the fol- 
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This range variables not one often encountered nature for the flat 
smooth bed regime. fact, probably not encountered often even labo- 
ratory flumes. nature the flat smooth bed regime seldom occurs depths 
less foot, and then only for velocities theorder fpsor great- 
er. Thus the should the order Actually from field ex- 
perience, the smallest known the writers the order 105. 

Perhaps the major appeal the Blasius formula simple power 
function, and easily handled mathematically. Thus, simplicity seems the 
mother invention. Other empirical formulas* have been offered well 
the Blasius describe the relation just beyond the critical Blasius 
considered sufficiently small range that first order approximation 
could used. Since actuality the relation not simple power function, 
the range can increased only using more complex equation. Her- 
man proposed, 


This equation would more difficult use than the inadequate 
that covers only the lower limit field data. Thus, either formula 
best only approximation, and perhaps shouldnot used for determining the 
powers and for flow open channels. 

The discharge formula proposed Mr. Liu was used compute the ve- 
locity for three reaches stream channel where measurements velocity, 
energy slope, and depths flow were available. comparisonof the observed 
and computed velocities the dune and smooth flat bed regimes shown 
Fig. 1(a). The measured velocity for all measurements higher the dune 
regime and lower the flat smooth bed regime than the velocity predicted 
the Liu equation. The range scatter the order 500%. 


Hydr. Engr., Geol. Survey, Washington, 
Hydr. Engr., Geol. Survey, Washington, 
Fluid Dynamics,” Prandtl, Hafner Pub. Co., New York, 1952. 
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(a) SLOPES MEASURED 
THE FIELD 
River 
| 
(6) SLOPES MEASURED FROM 
TOPOGRAPHIC MAPS 
Pigeon Roost Creek near Miss. 
Fivemile Creek near Wyo. 
.25mm. 
Ne 
MEASURED VELOCITY ft/sec. 


FIG, 1.—COMPARISON MEASURED VELOCITY WITH VELOCITY 
COMPUTED FROM LIU FORMULA 
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DISCUSSION 


similar plot for three additional reaches shown Fig. 1(b). For these 
reaches, the computed all flat smooth bedregime measurements 
are too high, and those for all but one dune regime measurement are too low. 
The slopes usedin computing these velocities fromthe Liu equation were taken 
from topographic maps. The predicted velocity appears better approxi- 
mation the measured velocity for the larger grain sizes. 

The criterion offered for prediction regime (Fig. paper) based 
upon V,/w This plot first was introduced closure paper 
Mr. and therefore has not received full discussion. Wd/v has replaced 
plot Liu previously offered. Mr. Liu states, “This permissible 


dune regime 
plane bed regime 


regime 


FIG, 2.—COMPARISON MEASUREMENTS KNOWN BED CONFIGURATION 
WITH LIU’S PREDICTED BED CONFIGURATION 


* 
because Actually this permissible only within the mean- 
ing dimensional analysis. Even there not entirely permissible, since 
and are not independent. mechanics upon which the criteria for vari- 
ous bed configurations are does not enter intoeither form except 
posteriori. Actually, there are six possible Reynolds numbers combining 
and with and The four involving and are identically acceptable 
empirically and equally unacceptable physical significance implied. The 
term chosen for ease use only. This field muddy one, and 


“Mechanics Sediment Ripple Formation,” Liu, Proceedings, ASCE, Vol. 
83, No. April, 1957. 


fi 
| ° + 
| | re) 
fe) 


mere manipulations dimensional analysis will not make any clearer. 
any rate, the criterion offered using V,/w and does not define the bed 
configuration natural channels. Data for the Rio Grande River near Ber- 
nalillo, New Mexico, and for Elkhorn River near Waterloo, Nebraska, are shown 
Fig. emphasize this point. 

Mr. Liu further states that the Froude number “can considered 
index indicating the energy loss caused the surface waves.” Since 


the implicationis that for given measuring site where the energy 
slope remains approximately constant, the energy loss caused surface waves 
function the Chezy For the flat smooth bed regime, Chezy 
the order twotimes that for the dune regime, and essentially equal that 
for the antidune regime. Therefore, follows that the energy loss caused 
the surface waves about for flat smooth bedand antidune regimes. 
The plane bed has surface waves, the antidune violent surface waves. Also, 
implied that the energy loss due waves for the flat smooth bed regime 
twice that for dune regime, although for the flat smooth bed regime both 
water surface and bed are believed the physical significance as- 
signed the Froude number may more apparent than real. 

The author should state whether the points used define Fig. are 
individual points, are averages series runs. From the small amount 
knowledge gained through working with datafor alluvial streams, does not 
seem that the data ever plot well indicated the figure. The scatter 
Figs. and herein seem more typical. 


GAMAL MOSTAFA.®—The authors have tackled important problem 
hydraulic engineering and have presented thorough study the present 
status research towards its solution. Attempts extend the applicability 
flow equations, which were originally derived for rigid boundaries, the 
field alluvial streams have far been unsuccessful. The authors have 
therefore turned empiricism and have proposed flow formula, Eq. 31, con- 
taining three sediment factors and which vary for different bed materi- 
and configuration. guide the application their method, they gave 
Fig. which can used when needed predict bed configuration. 

estimate and the authors gave Figs. and 12, which were de- 
duced from some collected data. The considerable jumps the values 
and when sand ripples appear dunes disappear may give the impression 
that suchchange bed configuration occur suddenly. this was true, 
thought that there should critical state for such sudden transfer. Con- 
sidering the movement sand bed under increasing flow traction, the 
change bed shape from plane ripple usually noticed just after the be- 
ginning general bed movement. The development sand ripples gradual 
and their transformation into dunes also gradual. This was first illustrated 
seventy years ago, who also noticed the gradual flattening the 


Dir., Research Sta., Delta Barrage, Egypt, V.A.R. 


“An Introduction Fluvial Hydraulics,” Leliavsky, Constable and Co., Ltd., 
London, 1954, 11. 
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DISCUSSION 


bed relatively high flow velocities. Straub,8 found that the maximum 
height dunes was reachedat critical flow conditions, corresponding mini- 
mum specific energy, and that the velocity flow rose above the critical, 
the undulations became less and less until finally the bed became smooth. 
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FIG, A.—NILE WATER PROFILES BETWEEN ASWAN AND ISNA 


Fig. gives the water surface profiles the Nile between the Aswan dam 
and Isna barrage six different discharges ranging from very high very 
low. During the field measurements, each case, the discharge was kept 


“Discussion Sand Movement Fluvial Models,” Straub, Transactions, 
ASCE, Vol. 61, 1935, 867. 
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constant for period long enough maintain steady flow. certain more 
less straight reach, long, located upstream the backwater effect 
Isna barrage, was chosen for testing the authors’ flow formula. Mechanical 
analysis bed material that reach gave mean sand diameter 0.25 mm. 
Table summary the results. the computations and which were 
deduced from Figs. and bed), were taken 0.44 and 0.31 respec- 
tively. The flood discharges corresponded the transition region Fig. 
and their computed were 11.61 and 11.40. For medium discharges was 
reduced 10.1 and the bed configuration was “dune.” low discharges 
rose again and reached value the lowest discharge. This was ex- 


plained the authors’ Fig. which showed that the bed condition was closer 
“ripple.” 


TABLE 1.—APPLICATION NILE DATA 


Discharge| Discharge 


water 
depth 


10900 
10200 
7500 
3360 
1620 
871 


Application the authors’ formula the Nile data has therefore proved 
satisfactory. Discharges estimated the formula differed from measured 
discharges more than 20%. Bearing mind that Manning’s formula, for 
example, may lead error over 100% the estimation river flow, 
the authors’ formula can certainly considered great improvement. 

would very interesting know how the proposed formula fits its 
application the Missouri River, which knownto carry heavy load sus- 
pension, and also its application alluvial streams with steep slopes. 


JAMES ASCE, and CARL NORDIN, au- 
thors are commended for analysis which covers wide scope the 
available flume and canal data, and for the presentation 
which, may adapted withease toeither design problems 
the prediction channel behavior. Moreover, the relationship between 
bed configuration and mean velocity flow appears related the phe- 
stage-discharge discontinuities which exhibit themselves true al- 
luvial channels, and the results the application the proposed formula 
data from measurements the Rio Grande New Mexico suggest that the 
formula may adapted some natural channel conditions. However, certain 
modifications are necessary due the fact that the authors’ Fig. does not 
adequately define bed configuration the field case. 

Current studies (1960) indicate that the true sand the Rio 
Grande exhibit stage-discharge discontinuities which apparently represent 

Hydr. Engr., Geol. Survey, Albuquerque, Mex. 

Hydr. Engr., Geol. Survey, Albuquerque, Mex. 


(1) (2) (3) (4) (5) (6) (7) (8) 
15.57 9.5 7.4 8.3 2.6 7.3 11.61 
14.57 9.3 7.0 7.95 11.40 
10.71 8.33 6.3 7.16 2.24 7.3 10.15 
4,80 5.8 5.27 1.64 7.3 
2.32 2.72 5.6 3.87 1.21 7.3 
1.24 5.1 2.66 .83 7.3 16.08 
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changes regimes flow, and which may detected plots hydraulic 
radius versus velocity. true alluvial sand bed channel used this 
discussion has the same characteristics the authors’ alluvial channel, that 
the bed material granular and cohesionless and the amount supply sedi- 
ment equal tothe amount sediment transport.) Fig. represents one such 
plot selected measurements the Rio Grande near Bernalillo, Mex. The 
flow may divided into lower and upper regime, with transition regime 


1.0 
0.5 


Velocity fps 


FIG, A.—RELATION HYDRAULIC RADIUS VELOCITY FOR RIO 
GRANDE NEAR BERNALILLO, NEW MEXICO 


between. The exact nature all changes which occur with 
changes regimes flow not known, but because the decreases resist- 
ance toflow which must accompany these changes, assumed that the lower 
regime represents flow over dunes and sand bars, and that during the tran- 
sition the upper regime, these bed features are obliterated flattened 
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the the flow. was found that the velocities for points which fall along 
the lower regime curve may computed using the authors’ curves for the sand 
bar configuration, and that velocities for points the transition stage between 
lower and upper regime may computed using the proposed curves for the 
flat bed configuration. Velocity measurements which plot along the upper re- 
gime curve Fig. appear represent flow conditions associated with flat 
bed antidunes and standing waves, and for this regime flow was possi- 
ble modify the curve actual field conditions, suggested the au- 
thors. 

total measurements was applied the proposed formula. These 
measurements covered range depthfrom about to9 ft, velocity from 


mm. Values and V,/w were computed using average median particle 


size and average water surface slope for the reach. These values were plot- 
ted the authors’ Fig. Eleven points fell the dune bed configuration and 
points fell the transition (flat) bed configuration. the points, only 
the values computed using the authors’ Ca, and were within 20% error. 
Classification into regime then proceeded the basis plots hydraulic 
radius versus velocity and field observations made the time the mea- 
surement water surface appearance andof relative bed compaction. Veloci- 
ties for these points were then computed using the authors’ values and 
and assuming that the lower regime points represent dune bed configuration 
and that regime points represent the flat About 
60% the values computed were within 20% error, showing marked im- 
provement over the classification bed configuration based upon V,/w and 

general, the field observations were value indicating where stand- 
ing waves were present and where the channel was hard and flat, which all 
cases were associated with the upper regime. The break the hydraulic 
radius-velocity curves was not well defined for all other stations for the 
Bernalillo station, and considerable difficulty was encountered differentiating 
the lower and the transition regime, especially where observations relative 
bed compaction and water surface appearance were not noted. 

For all other measurements which fell the upper regime, value 
was computed using the given exponents for and for plane bed and the 
measured velocities, and average values were used define the modified 
curve shown Fig. The lower portion the curve falls considerably be- 
low the values which might expected for these particle sizes, and was de- 
termined that for these stations portions the bed scoured down rock and 
gravel higher stages flow. Further, the writers believe that reaches 
the Rio Grande which have median particle sizes bed material greater than 
about 0.40 cannot true sand channel streams; that say, 
the bed configuration, velocity, and roughness parameters show unique re- 
lation with change regime flow. This fact limits the usefulness the 
curve those reaches which meet these criteriafor true sand channel streams. 

The results this study indicate that the proposed method for computing 
mean velocity flow can probably applied successfully true sand chan- 
nels such portions the Rio Grande, but that additional criteria beyond the 


relation V,/w are needed define adequately bed configurations. 
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The writers further any criteriafor determining regime flow 
true sand bed channels that based shear velocity cannot sufficient. 
energy gradient assumed for reach which has the hydraulic 
radius-velocity characteristics indicated Fig. then shear velocity be- 
comes function only depth flow hydraulic radius, and possible 
that more than one flow regime can exist for given hydraulic radius. 

conclusion, the writers propose that the following items included 
the authors’ suggestions for future research: 


The method should extended include the “antidune” bed configuration 
additional information becomes available. 

Additional field data are needed, with special attention directed the 
development equipment and techniques for measuring actual bed configu- 
rations. River profile studies such the one conducted Carey and 
shed considerable light upon the changes bed configurations 
which occur with change stage, and further investigations along this line 
would great value. 

The relation sediment discharge bed configuration, and its effect 
upon the proposed formula should investigated. 


Changes the Beds Alluvial Rivers,” Walter Carey and 
Dean Keller, Proceedings, ASCE, Vol. 83, No. August, 1957. 
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EFFECT AQUIFER TURBULENCE WELL 


effect turbulence head losses aquifers near 
pumped wells has been considered the author. the West many wells are 
constructed with envelope gravel prevent sand from moving into the 
well. the light this paper would interesting consider the effects 
turbulence flow gravel packs. 

the gravel pack both velocity and grain size are larger than the aqui- 
fer, therefore turbulent flow more likely exist. The permeability the 
gravel pack should several times higher than that the aquifer. The head 
loss resulting from turbulence, therefore, may not important. 

tests have been made the writer for the purpose determin- 
ing proper size and uniformity for gravel packs for use with wells various 
aquifers. The test model represented 30-in. diameter well with 12-in. 
diameter screen, size representative many irrigation wells. Seven pie- 
zometers along radius the gravel pack permitted head readings. Dis- 
charges ranged from the equivalent 100 gal per min per vertical 
well. Reynolds numbers were computed the method used the author. 
Reynolds numbers the gravel packs the well screen varied from 2.3 
94.0. Most were greater than 10. 

determine whether flow the gravel pack the model was lami- 
nar transition regime, the Reynolds number the well screen was calcu- 


lated. The value the expression was determined according 


the author’s table. Then using the head readings two points the gravel 
pack, one near the pack-aquifer interface and the other near the well screen, 
the value was calculated for two cases, listedunder the heading “Con- 
clusions from Analysis Test Data” and equal (assumed laminar flow). 
The head loss curve the gravel pack was plotted for each and 
compared with the experimental intermediate points the gravel pack. 

For Reynolds numbers the order 20, the assumed laminar flow curve 
agrees slightly better with the experimental data than does the transition curve. 
This illustrated Fig. for one test typical several this Reynolds 
number. The agreement the data with the laminar curve therefore ap- 
parently real rather than due experimental error. For tests with higher 
Reynolds numbers the well screen the experimental data generally plotted 
between the assumed laminar and transition curves. may concluded that 
some head loss due turbulence often exists gravel packs for Reynolds 


November, 1959, Joe Mogg. 
Engr., Western Soil and Water Mgt. Research Branch, Agric. Research Ser- 
vice, Dept. Agric., Fort Collins, Colo. 
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numbers greater than 20. The values velocity exponent listed the author 
are probably slightly high due the factor safety that the author has included. 

The head losses the model tests were in. less through properly se- 
lected and developed gravel packs. This value would include the component 
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FIG, 1.—HEAD LOSS CURVES FOR ASSUMED LAMINAR AND TRANSITION FLOWS 
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FIG, 2.—HEAD LOSS THROUGH WELL SCREENS 


loss due turbulence. For pumping lift more the total head loss 
through the gravel pack would amount less than the total lift and the 
component due turbulence could surely neglected. 

Reynolds numbers the aquifer the pack-aquifer interface were mostly 
less thanten and inall but two tests were less than 20. would seem that head 
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loss due turbulence the aquifer large diameter gravel-packed well 
can also neglected computing drawdown wells. 

Experimental data obtained byG. Corey, Jr.2 (Fig. support the author’s 
observation that screens with insufficient open area may cause significant head 
losses. Corey’s tests showed that head losses through the screen decreased 
the percentage open area, increased. For above some critical 
value, head losses remained constant. This critical value was dependent 
discharge and the size material surrounding the screen. Corey also found 
that the head losses through the gravel packs were quite small compared 
total losses through the well screens. 


“Hydraulic Properties Well Screens,” Gilbert Lee Corey, Jr., Unpublished 
Master’s Thesis, Colorado State Univ., June, 
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THE VORTEX CHAMBER AUTOMATIC FLOW-CONTROL 
Discussion John Stevens 


JOHN STEVENS, ASCE.—This paper was obviously submitted before 
Paper 1461 the author’s reference was printed Transactions Paper 
3004, Vol. 124, 1959, 871. Transactions, the Vortex Number the au- 
thors’ signs and symbols has been changed tothe Kolf Vortex Number and given 
the symbol recognition the original contribution the science hy- 
draulics Richard Kolf, ASCE. this discussion the symbol is, 
therefore, used instead 

The authors have added our knowledge regarding vortex flow through 
horizontal orifices. particular they have shown that the orifice discharge 
coefficient can determined approximately from the velocity the water 
entering the vortex chamber. They have also established the relation between 


and the angle divergence the water from the horizontal leaves 
the orifice, 


Otherwise had heretofore been determined from the circulation the orifice 
diameter and the average velocity water, through the orifice, 


stated the authors, the water surface profile computations were 
measurement made only one diameter away from the orifice” that is, 
substituting for Fig. was found the University Wisconsin ex- 
periments that coordinates this point were fair index the hyperbolic 
water surface profile free vortex2 given (28) except perhaps near the 
orifice where viscous forces are evidence. 

The writer was interested knowing how closely the the two 
methods agreed. From the authors’ table has been prepared. 


Cons. Hydr. Engr., Stevens and Thompson, and Leupold and Stevens Instruments, 


Inc., Portland, Ore. 
Transactions, ASCE, Vol. 124, 1959. 
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The average the ratios the last columnis 0.985 indicating that the 
age value either approach was within 1/2% each other, although 
the maximum variation was 11%. With more refined method measuring 
this angle, which the authors admit was difficult the Marquette University 
laboratory, better individual agreements could doubtless expected. 

comparison the discharge coefficients versus the Kolf vortex number 
for both the and the Wisconsin experiments are shown Fig. 


FIG, 1.—COMPARISON COEFFICIENTS AND VORTEX NUMBER 


Curve takenfrom the authors’ Fig. and defined experimental data, 
while Curve from the composite curve the authors’ reference the 
solid portion which similarly defined. The deflection toward Grandma’s 
near the bottom Curve might disappear some allowance could in- 


cluded for the non-hyperbolic portion the vortex surface profile water 
approaches the orifice. 
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THE SETTLING PROPERTIES 


Discussions Lucien Brush, Jr. and Hau-Wong Ho, and Hall 


LUCIEN BRUSH, and HAU-WONG HO.2—The author deserves con- 
siderable praise for his lucid presentation interesting subject. The de- 
tailed presentation the continuity equation for sediment-water suspension 
The experiment and the interpretations the results are clearly presented, 
are the comments similarity, and the speculations the application 
the multiple-depth analysis sediment problems involving open-channel 
flow and sediment removal. addition the interpretations made the au- 
thor concerning the experiment, several other interpretations exist which the 
writers feel are worthy comment. Because the core the author’s paper 
based the so-called measurement the local mean settling velocity, 
this discussion will directed toward the experimental and interpretive as- 
pects the investigation. 

presenting the continuity equation for sediment suspensions, the author 
states that turbulence will appear explicitly diffusion terms, where- 
flocculation hindered settling due turbulence will occur implicitly 
the fall-velocity term. Nevertheless, the results obtained the author for 
the multiple-depth analysis are interpreted due flocculation alone. 
Thus, any effect due turbulence, either implicitly explicitly, discount- 
ed. There seems little justification ignoring turbulence, because 
most certainly present, least initially, the procedure outlined the author 
followed making the analysis. Turbulence enters the system result 
“thoroughly mixing” the sample (to obtain uniform re- 
sult pouring the sample into the measuring cylinder. Obviously, the intensity 
the disturbance caused this experimental procedure will decrease with 
time due the viscous effects the fluid, but the net outcome the intro- 
duction decaying turbulence will generate curves similar those ob- 
tained the author Fig. 8—with without the presence flocculation. 

Although difficult measure the scale and intensity the introduced 
turbulence without running rather elaborate experiments, possible make 
qualitative estimate the amountof time that effect would 
importance experiment similar that made the author. The writers 
placed fluid droplets with specific gravity slightly greater than one con- 
tainer filledwith water and thoroughly mixed the liquid obtain dispersed 
suspension. The contents this container were then quickly poured into 
cylinder which had dimensions similar the one described the author. Ob- 
servations were then made the noting the direction 


December, 1959, Ronald Laughlin, Jr. 


Research Engr., Iowa Inst. Hydr. Research, and Asst. Prof., Dept. Mechanics 
and Hydr.; State Univ. Iowa City, Iowa. 
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motion the fluid droplets. Three minutes after pouring, was noted that 
about 30% 40% the fluid droplets were moving upward against gravity. 
After min, less than the particles were moving upward, but the effect 
was still noticeable. Thus, there seems reasonable evidence that some 
turbulence effect would exist for about min after pouring the 


would follow therefore, that some doubt might raised the applicability 
the author’s Eq. 


computing the local mean settling velocity, because diffusion term should 
included the equation. The appropriate equation would be, 


Integration this equation with respect will yield, 


instead Eq. 18, given the author. numerical solution the integral 
may obtained determining the slope, for particular depth, plot 
the area under the concentration curve against time, but the numerical solution 
will not discriminate between the author’s Eq. and Eq. Thus, the nu- 
merical result divided the concentration, apparent fall velocity—here 
called W,—will determined which will approach the author’s only the 
absence turbulence the turbulence decays sufficiently cause the 
diffusion term insignificant. Assurance that the diffusion term Eq. 
important initially, may obtained considering that order obtain 
initially uniform dispersion sediment, the following equation must apply 


wherein definition the two terms must equal. the turbulence decays, 
the diffusion term will become less important and will cause rise 
peak and then decrease function time much the same manner re- 
ported the author (Fig. for flocculation. follows therefore, that the ef- 
fect decaying turbulence computations the local mean settling velocity 
cannot readily distinguished from the effect caused flocculation without 
some other independent check. also note that the peak the 
curve shown the author (Fig. 8)occurs about min min after the in- 
itial pouring which the same order magnitude the time necessary for 
the turbulence effect become negligible, estimated the writers. 

Aside from the explicit effect turbulence, implicit effect may also oc- 
cur; that is, the settling individual grams will hindered—at decreasing 
pouring—as resultof turbulence. Decaying, hindered settling would 
also tend cause peak occur fall velocity time graph. order 
visualize these various effects, schematic illustration, showing plots fall 
velocity versus time for particular sampling depth given Fig. 
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Each the reported turbulence effects—decaying turbulence decaying 
hindrance settling—would tend cause the lines constant Z/t slope 
towards the axis concentration depth graph (Fig. and therefore be- 
have similar suspension which flocculation occurs. Because this and 
the previously stated arguments, there considerable doubt the minds 


settling with 


hindered settling 


Time 


FIG, TURBULENCE, HINDERED SETTLING, AND FLOCCULATION 
COMPUTATIONS THE LOCAL MEAN SETTLING VELOCITY. 


the writers whether the author’s experiment isolates the effect floccu- 
lation. Presumably, great depths, turbulence would decay before appre- 
ciable change the concentration gradient would occur, thereby reducing the 
importance including diffusion term the continuity equation. Perhaps 
this would enable the separation flocculation from the other effects, but 


order establish the proper depths for sampling, further experimentation 
would required. 


serves narrow the gap theoretical approach the 
process and practical design. 

common all the theoretical approaches from the time Hazen that 
basic assumptions the hydraulic flow tank and the settling character- 
istics the suspension have made which are necessarily simplified 
the point that practical application the resultant theory doubtful value. 

The settling characteristics the suspension, particular, factor 
which has not received the attention whichit warrants and the 
tions made that the particles have uniform settling rate, that they settle 
with uniform acceleration, far removed from the truth, particularly where 
sewage concerned. The graphs this paper clearly illustrate this point. 


Chf. Planning Engr., City Engr.’s Dept., Johannesburg, Africa. 
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The multiple depth pipette analysis suggested, deals adequately with this 
factor, but its application likely difficult the least laborious 
practice. 

the search for simple practical method determining the settling 
characteristics sewage, the writer has made close study the process 
takes place the laboratory and brief report the developments may 
interest. Ooservations show that even after sewage sample has been 
standing for hour two, flocs can seen settling such high rates that 
they settle from near the surface the bottom few minutes. 


presumed from these observations that settling sewage proceeds 
follows: 


Coarse material settles out with few minutes quiescence only. 

The fine material remains suspension slowly forming flocs. 

flocs reach critical size they commence settle, and their 
removal shallow vessel comparatively rapid. 


this basis, the selection experimental method determining the 
settling characteristics sewage, can made the assumptions that: 


The removal the coarse material sewage easily achieved and 
may ignored tank design (except far sludge storage and removal 
concerned). 

The fine material remains appreciable period while 
flocculation proceeds. 

Once flocs critical size are produced, their removal, compared 
the times required for their formation, instantaneous. 


Based these assumptions, therefore, samples taken from the surface 
quiescent sewage intervals settlement proceeds will reflect the settling 
characteristics that sewage. the application the method practice, 
may added that far the instantaneous removal flocs, once formed, 
concerned, the assumption can justified the overflow rate the tank 
under consideration less than the settling velocity the floc. 

The apparatus used the writer determine the settling characteristics 
sewage, consists drum in. diameter and in. deep which 
small floating draw off, 1/4 in. bore tube suspended two corks, fitted. 
batch sewage added and samples removed from the surface inter- 
vals. The samples are analysed either for suspended solids preferably for 
O.A.) this more significance works design. The experi- 
repeated many the results averaged and represented graphi- 
cally quality effluent (or surface liquor) against time. The method 
analysis well adapted demonstrate the value mechanical flocculation, 
the batch after adding test vessel may agitated the vessel for 
given time, given paddle speeds, after which the paddle mechanism re- 
moved, the floating draw off replaced and samples taken the usual way. The 
use two identical vessels simultaneously enables direct comparison 
made between pure quiescence and flocculation followed quiescence. Fig. 
shows typical settling characteristics analysis curve domestic sewage 
and the characteristics the same sewage when agitated oscillating 
paddle for min maximum paddle velocity 4.5 ips. 

would appear, first sight, easy determine the efficiency sedi- 
mentation tank relating its effluent quality time quiescence (or effec- 
tive detention period) reference the settling characteristics curve. Vari- 
ations settling characteristics and sewage strength during the day are such, 
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however, that reliable measure can only obtained averaging large num- 
bers readings. further factor that flocculation does not proceed 
fixed rate but greatly speeded slow mechanical agitation. Consequent- 
ly, must occur varying degrees all sedimentation tanks. For graph 
which direct reference could made, would, therefore, necessary 
either: 


Determine the extent and duration the agitation occurring the tank 
and duplicate this the vessel used for analysis 

optimum conditions for flocculation the vessel and thus 
obtain curve which could assumed represent the optimum and ideal con- 
dition. 


Time FLOCCULATION 


SPEED 


SAMPLE 


PV. (P. 


Time 


FIG, 1.—SETTLING CHARACTERISTIC ANALYSIS CURVE 


obvious that even with the simplified approach outlined previously, the 
prediction tank performance from basic data bristles with difficulties 
which unknown factors abound. hoped that further research such 
that this paper will carriedout tobring more light this fasci- 
nating problem. 
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TIDAL CHARACTERISTICS FROM HARMONIC 


Discussion Remig Papp 


REMIG paper fills gap our information about 
tides and the author should congratulated for it. 

would like add supplementary comment the section: Meteorological 
Effects. addition the winds, the barometric pressure can also influence 
the height the tide. After high barometric pressure the tide will lower 
than computed and predicted; after low pressure period will high- 
er. The difference could quite appreciable and may amount in. in. 

This writer has participated the port construction the Port Dun- 
querke, France, 1931-32. For the construction the Jetee Est, rein- 
forced concrete caissons, precast the inside the harbor, had towed 
out the end the jetty construction. Due the draft the caissons this 
towing could made only twice month during highest tide and every major 
deviation from the predicted tide level could have been critical. 

The field office has had registering barometer and the differences between 
the predicted andactual tide elevations were plotted under the barometric pres- 
sure curve. plotting the plus deviation downward and the minus deviation 
upward, the two curves became similar. The tide curve, however, had hada 
pronounced lag behind the barometric curve which might amount 10, 12, 
more hours. difficulty that the differences deviations the tide could 
observed the high and low points only other points are generally 
not computed, having importance. 

The notes the observations were lost during the war but they were not 
long enough draw too many conclusions from them. would be, however, 
quite easy establish this type observation few harbors with high 
tide range and evaluate the the barometric pressure tide ele- 
vation. conceivable that after collecting enough data, relation between 
the angle slope the barometric pressure curve and the deviation tide 
elevation can established for each local condition. 

narrow estuaries the wind effect and the influence the barometric pres- 
sure may have the same sign and may have superimposed they may 
have opposite signs. such cases the direction and force the wind should 
also taken into consideration. 


December, 1959, Bernard Zetler. 
Hazen and Sawyer Engrs., New York. 
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EARLY HISTORY HYDROMETRY THE UNITED 
Discussions William Hoyt and Shaw Cole 


WILLIAM Kolupaila uses the beginning World 
War (1914) the end the pioneer stage the field water measurement 
known stream gaging. order complete his listing pioneers this 
field seems appropriate include other engineers who added the know- 
ledge the field those early days. The writer refers the following who 
were district engineers the United States Geological Survey, Dept. the In- 
terior (USGS), 1918 and who were not otherwise mentioned the paper: 

Hartwell, New England and New Jersey; Horton, the Ohio 
River Basin and who initiated the first nationwide system collecting esti- 
mates developed and potential water power; William Lamb, Montana; 
McGlashan, California; Parker, Washington state and later chief 
hydraulic engineer (1939 1946); Paulsen, the south Atlantic states 
and later chief hydraulic engineer (1946 1957); Purton the Great 
Basin; Rice, Kansas; and Stevens, the Middle Atlantic region. 

Mr. Kolupaila suggests that those who measure water might well called 
“hydrometrists.” Perhaps the artof stream gaginghas progressed tothe point 
where such designation would appropriate. The engineers referred 
the paper listed herein were far more than hydrometrists. Had they been 
only what the definition implies, the Water Resources Division the USGS 
would not the authoritative agency national water resources today. 

Prior World War the actual time spent measuring streams took but 
each hr, fact towhich the writer can personally attestas one 
the pioneers listed Mr. Kolupaila. When not traveling train, livery, 
horseback, stream gagers were carpenters building structures; laborers 
chopping ice; public relation experts keeping farmers, ranchers, and others 
happy and reading gages twice day for maximum per month; comput- 
ing gage heights, preparing rating curves; making monthly reports Washing- 
ton; and above all, exercising expert salesmanship convincing state and city 
officials, heads irrigation districts and hydropower companies, and any 
others who would listen, the need for streamflow records. Without the fi- 
nancial support thus obtained and the enthusiasm engendered, the work would 
not have been maintained and certainly would not have continued ever- 
expanding rate the present time. With but small staffs and large regions 
cover, those charge the few district offices managed maintain some 
1,000 gaging stations total cost approximately $300,000, which 
one-half was provided state, public, and private agencies. 

some quarters these early engineers have been criticized for not taking 
time make myriad other observations deemed today essential for 


January, 1960, Steponas Kolupaila. 
Cons., Senate Select Committee Natl. Water Resources, Washington, 
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proper water management and also for not thinking more profoundly about the 
hydrology the river basins and the hydraulics and morphology the streams. 
general, can said that they were largely stream gagers, pure and sim- 
ple, but the records they recollected sticking their respective assignments 
were the beginning storehouse streamflow information that now covers 
some 200,000 station-yr records, with annual expenditure the USGS for 
all phases water information now running about $25,000,000, which co- 
operating bodies—local, state, and federal—still contribute about one-half. 

Unfortunately, the beginning the century there was immediate need 
(and the future need was not adequately recognized) for information about ground 
water, the quality water, the movement sediment and other hydrologic data 
now considered essential the management water resources. Even ex- 
tensive all the information water resources today, still far from 
sufficient meet present-day needs information now being collected the 
Senate Select Committee National Water Resources, which Senator Robert 
Kerr Oklahoma chairman, all indicative. 

The writer having opportunity analyze the extensive reports being 
received from the governors the states response request from 
Senator Kerr for information present water problems and those antici- 
pated between now and 1980, together with recommendations what should 
done toward their solution. Also, being reviewed the testimony some 
750 witnesses public hearings held throughout the country and numerous 
reports from federal agencies, consultants, and others. 

Although the writer somewhat familiar with the availability hydrologic 
data and deficiencies therein,2 has been truly amazed with the apparent 
need disclosed these sources for more hydrologic and meteorologic infor- 
mationof all kinds solutions the water problems today and those antici- 
pated 1980 are met. The extent which this information collected, 
analyzed, interpreted, and disseminated will have profound effect the ra- 


pidity and efficiency with which present and anticipated water problems will 
solved. 


pitot tube used measure the flow pipes. number published papers 
the subject have sought establish that there are two types, the simple and 
the combined. The simple pitot tube essentially bent tube with the nozzle 
orifice facing the flow. The total pressure against the orifice combi- 
nation velocity head plus static headand measure velocity pipe under 
pressure the static pressure must measured wall piezometer. This 
the type tube used Mills and Freeman. 

The combined type pitot tube included second tube that measured the 
static pressure. This was generally long tube that was suited open chan- 
nel flow measurements. The Pitometer was special type combined tube 
that used two identical orifices 180° apart, one facing upstream and the other 
downstream. This was developed and patented Edward Cole and giventhe 
coined name that Henri Pitot’s name. Its special features were its 
compact size that permitted installed easily pipe under pressure, 
thus making portable, its accuracy andits reversibility which provided 
simple means for checking zero. 


«water Facts For the Nation’s Future,” Langbein and Hoyt, Ronald Press Co., 
New York, 
Pres., The Pitometer Assocs., New York, 
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DISCUSSION 119 


Test Code for are important references for any one interested the 
development use the pitometer. 


«Pitot Tube Practice,” Edward Cole, Transactions, ASME, Vol. 57, 1935, 
281 
«Pitot Tubes Larger Pipes,” Edward Cole and Shaw Cole, Transactions, 
ASME, Vol. 61, 1939, 465. 
“Investigation Errors Pitot Tubes,” Hubbard, Transactions, ASME, 
August, 1939. 
“Power Test Code for Hydraulic Prime Movers,” 1949, pp. 135-165. 
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Journal the Hydraulics Division 


Proceedings the American Society Civil Engineers 


February, 1960 
41. The letter the end Eq. should made superscript. 


41. near the bottom the page the less than sign should 
moved closer the 1/2 the first line. 


47. the equation following Eq. the term should changed 


51. line the words “and the stage the floodplain” should 
changed “and the state the flood plain.” 
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DIVISION ACTIVITIES 
HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


NEWS 


JUNE 1960 


THE PURPOSE THE HYDRAULICS DIVISION 


“The advancement and dissemination scientific and engineering 
knowledge all branches hydraulics, hydrology, hydraulic engineer- 
ing and water resources. particular this shall embrace meteorology 
and hydrology the sciences dealing with the occurrence water 
the atmosphere, the earth surface and the ground, fluid mechanics 
for the understanding all flow phenomena, applied hydraulics for the 
design and planning hydraulic structures and comprehensive sys- 
tems, and those social, economic and administrative aspects basic 
the conservation and utilization water essential natural re- 
source.” 


All committees and task forces together with their members are listed 
the Official Register for 1960. 


Committee Meetings 


Division Technical Committees are expected extremely active the 
future. Several meetings have been held and others are being planned. Mem- 
bers the Hydraulics Division can look forward many interesting programs 
result this activity. Brief notes the committee meetings will ac- 
quaint you with some the activities. 


Committee Sedimentation 


This committee met Chicago April and 10, with Chairman 
Banks attending. The committee firmed the program for the session the 
Seattle Conference (see Seattle program this issue the Newsletter), 
made plans complete the Sedimentation Manual early date, laid plans 
for sessions future conventions and reviewed the activities existing task 
forces. Chairman Carlson will plan session for the April 10-14, 1961 


Note.—No. 1960-19 part the copyrighted Journal the Hydraulics Division, 
Proceedings the American Society Civil Engineers, Vol. 86, No. June, 1960. 


Copyright 1960 the American Society Civil Engineers. 
1960-19--1 
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Phoenix Convention under the theme “Scour Cohesive Materials,” and mem- 
ber Eagleson will plan session for the August 16-18, 1961 Urbana Hydraulics 
Conference under the theme “Reservoir Sedimentation.” The committee’s 
future plans include setting task force prepare bibliography 
translations and foreign literature various aspects sedimentation. 


Committee Hydraulic Structures 


This committee met March during the New Orleans Convention with 
and Hands attending. The activities the Committee and the eight task 
forces were reviewed detail and plans made for activities including ses- 
sions for future conventions. The aims, accomplishments, status, objectives, 
progress, and action for each task force were considered individually and 
detail. Extensive future activity expected. Many the task forces will 
plan and prepare material for future sessions. The program for the Seattle 
Conference has been completed (see Seattle program this issue the 
Newsletter). session for the Boston Convention (October 10-14) under the 
theme “Hydraulics Bridges” being arranged Hands. session 
for the Phoenix Convention, April 10-14, 1961. Sessions with themes 
“Hydraulics Culverts,” “Energy Dissipators,” and “Hydraulics Bridges,” 
are planned for the Urbana Hydraulics Conference, New York Conven- 
tion, and Houston Convention. The committee proposes set new task 
force “Hydraulics Navigation Locks” gather and present information 
recent progress and changes design practices. 


Committees Hydrometeorology, Surface Water Hydrology, and Flood 
Control 


These committees held joint meeting April and Chicago, with 
and Waananen the Hydrometeorology Committee, Chairman 
Foster, and members Ackermann, McDonald and Jens the 
Surface Water Hydrology Committee, Chairman Lenz the Flood Con- 
trol Committee, and Riesbol, Water Resources Coordinator attending. 
The committees planned and coordinated their program work for the year. 
Objectives and scope work were agreed upon and task assignments were 
made committee members. Sessions for the Society Convention Phoenix, 
Arizona, April 10-14, 1961, and the Hydraulics Division Conference 
Urbana, Illinois, August 16-18, 1961, were planned. Task group objectives 
and assignments were reviewed. 


Hydrologic Materials Preserved 


Because the increased interest and activities the Hydraulics Division 
Water Resources, recent article from the February 1960 issue “West- 
ern Water News” quoted below. The new Water Resources Center Archives 
unique the United States and should interest and value ASCE 
members interested water resources. 
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“ARCHIVES PROGRAM PURSUES HYDROLOGIC MATERIALS 


Center Seeks Cooperation Water Agencies 


part state-wide University California research program 
water resources, the Water Resources Center Archives the Berke- 
ley Campus gathering materials and data the water resources 
field. The collection will include all aspects water: water natu- 
ral resource and its utilization; irrigation; flood control; municipal and 
industrial water uses and problems; water rights; and water develop- 
ment projects. The emphasis material relating the State 
California. 

Besides actually acquiring material, the Archives interested 
locating and making biographical listings important collections which 
remain public private agencies and engineering offices. com- 
prehensive catalog will thus developed for the use researchers, 
describing and giving the location existing unique hydrologic data 
the State. 


Cooperation Sought 


Because the materials found the most valuable are not available 
the sense that they can purchased such most libraries purchase 
books, the Archives relies upon the cooperation interested individuals 


and agencies for augmenting its collection. Engineers’ reports and 
studies (unpublished materials, manuscripts, etc.), agency reports, 
pamphlets, hearings, legislation, specifications, maps, scrapbooks, 
newspaper clippings, speeches, campaign material, and books pertain- 


ing water: all are being sought. Interests include both current mate- 


rials and the historical aspects. 
The entire collections several men prominent the water re- 


sources engineering fields are now housed the Archives: among these 
important gift collections are those Frank Adams, Bernard 
Etcheverry, Charles Gilman Hyde and Baldwin Woods. 

Seven reports have been prepared and published the Archives, all 
them available the interested public: 


The Etcheverry Collection Water Resources Archives 
(June, 1958). 

Theses Engineering, Economic, Social and Legal Aspects 
Water (October, 1958). 

Watershed Management Research Data, Department 
Agriculture, Forest Service, California Forest and Range Experiment 
Station, Berkeley and Glendora, California (February, 1959). 

Publications and Reports Charles Gilman Hyde (July, 1959). 

Water Pollution Data, Regional Water Pollution Control Boards, 
State California (July, 1959). 

Bachelor Science Theses Water Resources Engineering, 
University California, Berkeley (August, 1959). 

Theses Water Resources, Stanford University, California 
Institute Technology, and University Southern California (Au- 
gust, 1959). 
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Agencies dealing with water any its forms uses are request- 
add the Archives their mailing lists. Persons having materials 
which have served their purposes are urged consider the Archives 


repository for the filing these materials. 

The Archives serves researchers the University community 
well the interested public. The activity under the direction 
Professor Johnson, Director Hydraulic Laboratories, Univer- 
sity California, Berkeley. Gerald Giefer the Librarian directly 
charge the Archives. 

The mailing address the Archives is: Water Resources Center 
Archives, Room Mechanics Building, University California, 
Berkeley California.” 


NINTH HYDRAULICS DIVISION CONFERENCE, SEATTLE, WASHINGTON 
August 17-19, 1960 


Include this conference your vacation plans. The University Wash- 
ington and the Seattle Section are cohosts and the setting will the Univer- 
sity Washington campus. 


The following program indicative the many interesting activities be- 


ing planned. 


General Information 


Housing and Dining Facilities Dormitory accommodations will available 


the University. Many hotels and motels are located within five mile 


radius the campus, and information concerning these will sent from the 


local planning committee upon request. University dining facilities will 
available for all meals. 


Social Events 


Pre-Conference Skagit River Tour, August This guided tour through 


the Seattle City Light integrated three-dam hydroelectric project, located 
the northern Cascade Mountains near the Canadian border. The all-day trip 
combines technical interests with rugged mountain scenery, and includes 
rides bus, boat, and inclined railway and lunch served the project. 


Salmon Barbecue, August Like the Skagit trip, this also planned 


family event. Scheduled for the evening the first day the formal con- 
ference, provides time for informal socializing boat trip through 


the Chittenden locks and across salt Puget Sound, followed typical 


regional salmon barbecue water. 


Evening Banquet, August This traditional event will the final formal 
social event the conference. 


Ladies’ and Children’s Programs program including teas, tours, and 


shopping being planned for the ladies. Activities for younger members 
the family are also being set up. 
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Technical Program 


Session Hydraulic Structures: 


Presiding: Dickinson, Vice Chairman, Executive Committee and 
Peterka, Chairman, Hydraulic Structures Committee. 


Effects Cavitation and Negative Pressures the Passage 
Fish Hydraulic Turbines 

Glen Von Gunten, Chief, Planning and Reports Branch, Corps 
Engineers, Walla Walla, Washington 


Problems Associated with Water Temperature and Fish Passage 
Milo Bell, Professor, College Fisheries, University 
Washington, Seattle, Washington 


Fish Protective Facilities Tracy Pumping Plant, Central 
Valley Project, California 

Rhone, Bureau Reclamation, Hydraulic Laboratory 
Branch, Denver, Colorado; Daniel Bates, Fishery Research 
Biologist, Bureau Commercial Fisheries, Portland, Oregon. 


Chief Joseph Dam Stilling Basin—Damage and Analysis Con- 
tributing Factors 


Robert Gedney, Chief, Planning Section, Corps Engineers, 
Seattle, Washington. 


Session Tidal Hydraulics: 


Presiding: Eugene Fortson, Jr., Member, Executive Committee and 
Haywood Dewey, Jr., Chairman, Tidal Hydraulics Committee. 


Prototype Measurements Columbia River Estuary 
John Lockett, Army Engineer Division, North Pacific; 
Harold Kidby, Army Engineer District, Portland, Oregon. 


Detailed Current Studies the Puget Sound Region 

John Dermody and Dr. Maurice Rattray, Department Ocean- 
ography, University Washington, Seattle, Washington. 

Calibration Method for Hydraulic Model Subjected Two 
Tidal Entrances 


Board Canada, Pacific Oceanographic Group, Nanimo, 


Session Hydrology: 


Presiding: Riesbol, Water Resources Coordinator, and Franklin 
Snyder, Member, Hydrology Committee. 


Radar Observations Topographic Effects Precipitation 
the Pacific Northwest 


Herbert Kershaw, Jr., Wilfred Jensby, and Professor Fred 


Decker, Science Research Institute, Oregon State College, 
Corvallis, Oregon. 


Errors Measurement Reservoir Water Levels 


Walter Parsons, Jr., Corps Engineers, Sacramento, 
California 
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Fluctuations Annual River Flows 


Dr. Vujica Yedvjevich, Geological Survey, Washington, 


Effect Time and Space Distribution Rainfall the Runoff 


Hydrograph 
Amorocho and Professor Orlob, Department Civil 


Engineering, University California, Berkeley, California. 


Session Flood Control: 


Presiding: 


Session 


Presiding: 


Session 


Presiding: 


Kindsvater, Past Chairman, Executive Committee, and 
Francis Christian, Member, Flood Control Committee. 


Development Spillway Design Flood for Bruces Eddy Dam, 
North Fork Clearwater River, Idaho 


Melvin Ord and Robert Conway, Corps Engineers, Walla 
Walla District, Walla Walla, Washington. 


Seasonal Runoff Forecasting the Columbia River the 
Coastal Flow Index Method 

David Rockwood, Corps Engineers, North Pacific Division, 
Portland, Oregon. 


Flood Control Economics, Lower Columbia River 


Kenneth Case, Corps Engineers, Portland District, Portland, 
Oregon. 


The Hydrology Howard Hanson Dam and Reservoir 
Norman MacDonald, Corps Engineers, Seattle District, 
Seattle, Washington. 


Hydromechanics: 


Dr. Ippen, Chairman, Executive Committee, and Dr. 
Norman Brooks, Member, Hydromechanics Committee. 


Report Task Force Aearted Flow Open Channels 

Frank Campbell (Chairman Task Force), Chief, Hydraulic 
Analysis Branch, Waterways Experiment Station, Vicksburg, 
Mississippi. 


Symposium Mechanics Ground Water Flow 
Drainage Wells Two-Layered Aquifer 
Jacob, Groundwater Consultant, Los Angeles, California 


Nonsteady Flow Toward Wells Partially Penetrating Water- 
Table Aquifer 

Mahdi Hantush, New Mexico Institute Mining and Technology, 
Socorro, New Mexico. 


Salt Intrusion into Coastal Aquifers 
Professor Harold Henry, Michigan State University, East 
Lansing, Michigan 


Sedimentation: 


Martin, Secretary, Executive Committee and Paul Benedict, 
Member, Sedimentation Committee. 
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Formulas for Total Sediment Discharge Alluvial Streams 
Dr. Vito Vanoni, California Institute Technology, 
Pasadena, California. 


Transportation Materials Fluids Pipe 

Professor Chamberlain; Yotsukura, Research Assis- 
tant; Karaki, Assistant Professor; Dr. Albertson, 
Colorado State University, Fort Collins, Colorado. 


Stream Control Planning and Structures 
Bondurant, Hydraulic Engineer, Corps Engineers, 
Missouri River Division, Omaha, Nebraska. 


Legal Aspects Sedimentation 
Busby, Water Program Specialist, Soil Conservation 
Service, Berkeley, California. 


Other Meetings 


Don’t Forget the ASCE Boston Convention, October 10-14, 1960, which will 
include sessions Tidal Hydraulics, Hydrology, Hydromechanics, Sedimen- 


tation, and Hydraulic Structures. More this meeting the August News- 
letter. 


Geological Society America Annual Meeting 
(This repeat from the April Newsletter) 


The Division Engineering Geology GSA will have symposium 
ground water Denver, Colorado, October and 31, 1960, along with 
special papers related subjects. field trip major projects and points 
interest engineering geologists the Rockies west Denver planned 
for Saturday, October 29. expected that the new Water Resources Com- 
mittees the Hydraulics Division will find many common interests with the 
Division Engineering Geology GSA. 


Officers for the Division Engineering Geology include three ASCE members 
and are: 
Chairman, Safford Happ, AEC, Grand Junction, Colorado 
Vice Chairman, George Kiersch (F. ASCE), Southern Pacific Company, 
Oakland, California 
Secretary, Thomas Thompson (Aff. ASCE), San Francisco, California 
Counsellor, George Schwartz, University Minnesota, Minneapolis, 
Minnesota 
Robert Legget (F. ASCE), National Research Council 
Canada 


For Your Calendar 


ASCE Meetings 


August 17-19, 1960 Hyrdraulics Division Conference, 
Seattle 


October 10-14, 1960 ASCE, Boston Convention 
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April 10-14, 1961 
August 16-18, 1961 


October 16-20, 1961 
February 1962 

May 1962 

October 15-19, 1962 


Non-ASCE Meetings 
September 1960 


October 29-November 1960 


June 26-July 1961 


June 18-21, 1962 


preciated. 


Deadline dates for Newsletter contributions: August 1960 issue—June 20; 


October 1960 issue—August 20. 


Use the Hydraulics Division Newsletter 


You are urged continue use the Division Newsletter for announce- 
ments, inquiries, personnel news, committee reports, surveys and other 
items interest Division members. short note summarizing the high- 
lights committee meetings particularly requested. Suggestions for im- 
provement the Newsletter will appreciated. All contributions are ap- 


James Ball, Editor 


June, 1960 


ASCE, Phoenix Convention 


Hydraulics Division Conference, 
Urbana, Illinois 


ASCE, New York Convention 
ASCE, Houston Convention 

ASCE, Omaha Convention 
ASCE, Detroit Convention 


Hydraulic Turbine Research Sym- 
posium, IAHR, Nice, France 


Annual Meeting Geological Society 
America, Division Engineering 
Geology, Denver, Colorado 


7th Congress International Committee 
Large Dams, Rome, Italy 


4th National Congress Applied 
Mechanics, University California, 
Berkeley, California 


Box 7416 
Lakewood 15, Colorado 
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